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Abstract 

The maximal center gauge, combined with center projection, is a 
means to associate Yang-Mills lattice gauge configurations with closed 
center vortex world-surfaces. This technique allows to study center 
vortex physics in lattice gauge experiments. In the present work, the 
continuum analogue of the maximal center gauge is constructed. This 
sheds new light on the meaning of the procedure on the lattice and 
leads to a sketch of an effective vortex theory in the continuum. Fur- 
thermore, the manner in which center vortex configurations generate 
the Pontryagin index is investigated. The Pontryagin index is built up 
from self-intersections of the vortex world-surfaces, where it is crucial 
that the surfaces be globally non-oriented. 
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1 Introduction 



The infrared sector of strong interaction physics is characterized by nonper- 
turbative phenomena such as confinement, the breaking of chiral symmetry, 
and the Ua(1) problem. These phenomena are believed to be described by 
quantum chromodynamics; however, in order to understand in more detail 
how they arise, it is helpful to isolate the particular degrees of freedom of 
the theory which induce these phenomena. Chiral symmetry breaking and 
the Ua(1) problem can e.g. be successfully explained in terms of instanton 
effects HJ, whereas confinement eludes a description via instantons [[J. 

Viable explanations of confinement can be based on magnetic degrees of 
freedom such as Abelian magnetic monopoles [[| or center vortices ]§,[§]; 
considerable evidence in favor of such scenaria has accumulated in recent 
years from numerical lattice studies. Magnetic monopoles emerge in Yang- 
Mills theory as gauge fixing defects in the so-called Abelian gauges 0, where 
the Cartan subgroup of the gauge group under consideration is left unfixed. 
In all Abelian gauges considered, the vacuum behaves as a dual superconduc- 
tor Jt|; in particular, lattice calculations performed in the maximal Abelian 
gauge 0.0 exhibit Abelian and monopole dominance for the string tension. 
On the other hand, lattice calculations performed in the so-called maximal 



center gauge ||10|| , which will be discussed in detail in this work, have pro- 
vided evidence in favor of the notion that confinement is induced by center 
vortices. The vortex content of gauge field configurations produces virtually 
the full string tension; conversely, the string tension disappears in the ab- 
sence of vortices This property of center dominance persists at finite 



temperatures [ll|.[12], and the deconfinement phase transition can be un- 
derstood as a transition from a percolating vortex phase to a phase in which 
the vortices cease to percolate. 

The aforementioned magnetic degrees of freedom inducing confinement can 
furthermore be connected with the topological properties commonly thought 
to be carried by instantons; they may therefore ultimately be suited to pro- 
vide a unified picture of the different nonperturbative strong interaction phe- 
nomena mentioned further above. The magnetic monopoles of the maximal 
Abelian gauge have been empirically found to be correlated with instantons 
in lattice calculations |13| . In another Abelian gauge, the Polyakov gauge, the 



Pontryagin index can be entirely reconstructed from the magnetic monopole 



content of gauge field configurations Jl4|]. On the other hand, one of the 



foci of the present work is to elucidate the manner in which center vortex 
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configurations generate the Pontryagin index; it will be shown that the latter 
is produced by self-intersections of vortex world-surfaces, where it is crucial 
that the surfaces be globally non-oriented. In a recent lattice study [[y| , the 
orientability of the vortex surfaces was studied, with the result that these sur- 
faces are non-orient able in the confining phase and have non-trivial genus. 
Furthermore, in an ensemble devoid of vortices, chiral symmetry breaking 
disappears and all configurations belong to the trivial topological sector []16 



For the abovementioned lattice studies of vortex properties to become possi- 
ble, techniques had to be developed which allow to associate a given lattice 
gauge configuration with a configuration of vortex world-surfaces. Such tech- 
niques have only become available in relatively recent times; this has revived 
interest in the vortex picture, which received rather little attention after 



some early efforts following its inception [j3]l , , [ l7| . The techniques referred 
to here are the maximal center gauge and center projection |rO||,||lS|| already 
mentioned further above. While not all questions concerning the numerical 
stability of this procedure and its ability to faithfully represent thick center 
fluxes present in full gauge configurations have been settled [jnj, it has fa- 
cilitated a number of observations concerning vortex physics in addition to 
those already mentioned further above. E.g., the density of the center vor- 
tices obtained via the maximal center gauge and center projection displays 
the proper scaling behavior required under the renormalization group and, 
consequently, these vortices can be considered as physical objects ]20j (note 



erratum in |TT[), see also [ZTIj. Also binary correlations of vortex intersection 



points found on a given space-time plane display such scaling p2[ . Further- 
more, center vortices can account for the Casimir scaling behavior of the 
adjoint string tension if one associates a finite physical thickness with the 



vortex world- surfaces |23],p3|. In an indirect version of the maximal center 
gauge, which proceeds via an intermediate maximal Abelian gauge fixing, 
both Abelian magnetic monopoles and center vortices were detected JTj| ; 
the monopoles were found to be located on center vortices like beads on a 
string. 

One of the foremost aims of this work is to develop a continuum formulation 
for the center gauge and center projection procedure. This will also shed some 
additional light on issues concerning the procedure on the lattice. It should 
be clear that the attempt to construct an analogue of the maximal center 
gauge in the continuum encounters some conceptual obstacles not found in 
continuum formulations of Abelian gauges, which are quite straightforward. 
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These obstacles are related to the fact that the continuum theory is con- 
structed in terms of a gauge potential, defined in the adjoint representation 
of the gauge group, as opposed to the lattice formulation in terms of group 
elements on lattice links. The continuum gauge fields are thus insensitive to 
the center of the gauge group. It is a priori not clear how to extend the suc- 
cessful concept of center projection to the continuum theory. Nevertheless, it 
will prove possible in this work to derive a continuum analogue of the lattice 
procedure. 

In detail, the plan of this paper is as follows. In section |2|, the maximal center 
gauge and the emergence of center vortices after center projection is carefully 
analyzed on the lattice. The effect of this gauge on sample gauge fields of the 
vortex type is studied empirically; subsequently, the gauge fixing procedure 
on the lattice is reformulated in a manner which explicitly separates the 
Z(N) center and SU(N)/Z(N) coset degrees of freedom. The continuum 
limit of center and coset lattice configurations is considered in section [| 
along with the gauge transformations which arise in the course of maximal 
center gauge fixing. Explicit continuum representations of center vortices are 
given. This allows the construction of the continuum analogue of the maximal 
center gauge, the properties of which are then discussed extensively, including 
alternative formulations which allow to extract more general objects of the 
vortex type; for reasons mentioned further above, vortex surfaces made up of 
patches of different orientation are of particular interest, as well as thickened 
vortices. In section £|, the implementation of the maximal center gauge in 
the continuum path integral is considered, and a sketch of an effective vortex 
theory is given. Finally, the Pontryagin index generated by center vortex 
configurations is studied in section The concluding section ^] discusses 
some open issues to be settled in subsequent work. Some mathematical 
details are presented in the appendices. 
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2 Maximal Center Gauge and Center Vor- 
tices 



2.1 The Maximal Center Gauge 

Maximal center gauge is a gauge fixing condition defined in the lattice for- 
mulation of SU(N) Yang-Mills theory. It prescribes that the gauge freedom 

U,(x) - U^(x) = V(x + e M )t^(x)y+(z) (1) 

be used to transform the link variables U^x) of a given lattice configuration 
as close as possible to center elements of the gauge group, in a sense still to be 
specified. In fact, there exist many different specifications, or realizations, 
of the general idea formulated above. One possible explicit gauge fixing 
condition, first introduced in |18 for the case of 577(2) color| as the "Direct 
Maximal Center Gauge", is 

max£|tr[/J(x)| 2 (2) 

X,fJ, 

for arbitrary local SU{2) gauge transformations V. 

In order to make this more transparent, represent the SU(2) link variables 

as 

U — Co + icr (3) 

where Tj denotes the Pauli matrices and the parameters Co, Cj, i = 1,2,3, fulfil 
the constraint 

eg + c 2 = 1 (4) 

to guarantee unitarity of the links. Note that (§]) defines a three-sphere in 
the space of parameters c , q. In the following, the set of group elements 
with Co > will sometimes be referred to as the "northern hemisphere" of 
the 577(2) group, and the set of group elements with Co < as the "southern 
hemisphere". The poles of the sphere, Co = ±1, correspond to the center 
Z{2) of the 577(2) group. Correspondingly, the gauge condition (Q) reduces 
to 

^EtiW) ( 5 ) 



5 Generalizations and applications to SU(3) color were discussed in |2lJ],[E5| 
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which can be straightforwardly visualized in terms of a bias of the link vari- 
ables towards the poles of the SU (2) three-sphere. 

An important property of the condition (0) is the fact that it leaves the 
center part of the gauge freedom, i.e. (for the general case of a SU (N) gauge 
group) the Z(N) part, unfixed. There is no bias distinguishing between 
different center elements due to the absolute value prescription in (0). This 
is analogous to the so-called Abelian gauges || , which leave the Abelian part 
of the gauge freedom, i.e. the [^(l)]^ -1 part, unfixed. 

The specific condition (0) is only one of many ways to bias link variables 
towards the center of the gauge group. In fact, any condition of the type 

m £ x 5^(| trf/ /T(>)|) ( 6 ) 

with a monotonously rising function g implements the general idea of maxi- 
mal center gauge fixing. 

After transforming a lattice gauge configuration to the maximal center gauge, 
one can carry out a second step, called center projection, which serves to 
define the Z(N) vortex content of the configuration. Center projection cor- 
responds to mapping every link variable U onto the center element closest to 
U. In the case of the £77(2) gauge group, this simply means mapping 

U -> sign tr U , (7) 

i.e. 

c -> sign c . (8) 

In other words, all link variables located on the northern hemisphere are 
mapped onto the north pole of the SU (2) group, all links on the southern 
hemisphere onto the south pole. For general SU(N) group, center projec- 
tion leaves one with a Z(N) lattice, which in the canonical fashion can be 
associated with (D — 2)-dimensional vortex surfaces (where D denotes the 
space-time dimension) on the dual lattice: If any plaquette on the Z(N) 
lattice evaluates to a nontrivial center element, it is considered to be pierced 
by a vortex carrying the appropriate (quantized) flux; if it takes the trivial 
value unity, no vortex pierces the plaquette. The vortices defined in this way 
represent (D — 2)-dimensional surfaces in D space-time dimensions. The 
center flux carried by the vortices is continuous; there are no sources or sinks 
for this flux. In the case of a Z(2) lattice, this corresponds to the vortex 
surfaces being closed. In the general Z(N) case, the surfaces may branch, 
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while respecting flux conservation, since more than one nontrivial center flux 
exists. However, precisely due to flux conservation, one may think also of 
such branched surfaces in the Z(N) case as being composed of simple closed 
surfaces which in general share whole surface patches in space-time. 

2.2 Center vortices on the lattice 

For the purpose of illustrating the effect of the center gauge condition, con- 
sider the following sample continuum SU(2) gauge field configurations on a 
two-dimensional plane. Let the plane be parametrized in polar coordinates 
r, ip and define the gauge field as[] 

a v = ±^T 3 , a r = (9) 

This configuration, and the ones to follow below, can be regarded as slices 
of configurations defined in more space-time dimensions, and which, for the 
present purposes, can be considered constant in the additional dimensions. 
Also the Lorentz components of the gauge field associated with those addi- 
tional directions should be considered to vanish. The above configuration 
is a particular example of what in the following will be called a thin vor- 
tex (positioned at the origin of the plane). A general representation of thin 
vortices will be given in eq. (§2|). The defining property of a thin vortex 
is that any Wilson loop which encircles the vortex, i.e. in the present case, 
the origin, takes a nontrivial value in the center of the gauge group. For the 
case of SU{2) color, this corresponds to the value (—1). On the other hand, 
any other Wilson loop, not encircling the vortex, takes the value (+1). In 
other words, the vortex contributes a nontrivial center element of the gauge 
group to any Wilson loop it pierces. Correspondingly, the field strength is 
singular and localized at the position of the vortex. Note that both choices 
of sign in (||) lead to these properties. Thus, one can distinguish two dif- 
ferent orientations of vortices, associated with the direction of (magnetic) 
flux they describe (flux which is dual to Wilson loops will sometimes be re- 
ferred to as "magnetic" in the following). On the other hand, also (odd) 
integral multiples of the configurations (|9]) satisfy the defining properties of 

4 Gauge fields in the following are decomposed using antisymmetric generators T a of 
the gauge group satisfying [T a ,Tj,] = f a bcT c , with f a b c denoting the structure coefficients, 
and the normalization TrT a Tb = —5 a b/2. For SU(2), T a — —ir a /2, where r OJ a = 1,2,3, 
denotes the Pauli matrices. 
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center vortices. However, such configurations will not be considered as dis- 
tinct types of vortices, but simply as superpositions of vortices of the smallest 
Abelian magnetic flux required to induce the center element in question in a 
Wilson loop. 

Consider furthermore the more general configuration 

fir) 

A v = ±^T 3 , A r = (10) 

with /(0) = and /(oo) = 1. This configuration represents a thick vortex; 
in order for a Wilson loop to take the value (—1), it must circumscribe the 
vortex at a sufficient distance, the thickness of the vortex, determined by the 
form of the profile function fir). Also the field strength or the flux associated 
with a thick vortex are spread out as compared to the singular localized field 
strength of a thin vortex (^|). For certain purposes, such as calculating field 
strengths, it may be necessary to specify more precisely what one means by 
the singularity of the thin vortex (|9[) at the origin; this becomes relevant e.g. 
in Appendix 0. In such the thin vortex should always be viewed as 



the limit in which f(r) — > 1 for r ^ 0, but still f(r = 0) = in (10). 



These configurations can be put on a two-dimensional lattice in the canonical 
fashion by associating the path-ordered exponential 

U^x) = Pexp (- J p Ap{z)dz^ (11) 

with a link U^x) describing a path P originating from the point x and run- 
ning in the //-direction. Note that path-ordering V is not actually necessary 
for the Abelian configurations ([|) and (TO). This procedure of putting the 
continuum configurations on the lattice preserves the values of Wilson loops 
taken in the corresponding configurations, as described above. 

In practice, it is desirable to treat a lattice of finite extent with periodic 
boundary conditions, i.e. a lattice torus. In order to fulfil such boundary 
conditions, it is necessary to have zero net flux piercing the plane under 
consideration^. For this reason, in the following, linear superpositions of 
pairs of (thin or thick) vortices of opposite orientation will be examined, 

5 Consider a Wilson loop running around the edge of the lattice. Due to the periodic 
boundary conditions and the Abelian character of the configurations treated here, contri- 
butions from opposite sides of the loop cancel and the Wilson loop must take the value 
(+!)■ 



cf. Fig. [I] (a). Due to the Abelian character of the configurations defined 
above, the phases contributed to Wilson loops by single vortices in a many- 
vortex superposition simply add, without interference from path-ordering 
effects. Furthermore, in calculating the lattice links from the gauge fields, 
also the contributions from periodic images of the vortices must be taken 
into account! 5 !. 
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Figure 1: Two thin SU(2) vortices, cf. eq. (Q), placed on a space-time lat- 
tice (a), where the arrows in the lower left-hand corner specify the direction 
of the line integrals defining horizontal and vertical links, respectively. The 
open and filled circles denote the locations of the vortices, the difference in 
the symbols signifying that the vortices have opposite orientation. In (a), all 
lattice links deviate from unity, but all plaquettes multiply to (+1), except 
the plaquettes pierced by vortices, which multiply to (-1). The links mak- 
ing up those latter plaquettes take the values U = exp(T 3 7r/2) and V = W 
(up to small corrections from the respectively other vortex). After maximal 
center gauge fixing (b), all lattice links take the value (+1), except the ones 
displayed in fat print, which take the value (-1). There are alternative gauge 
fixing images, such as (c), for which the center gauge fixing functional is 
degenerate, and which are related to the configuration (b) by gauge trans- 
formations purely from the center of the gauge group. 



6 In practice, for better convergence of this sum, the authors covered the plane with 
positively and negatively oriented vortices in checkerboard fashion. Therefore, the lattice 
torus actually used in numerical computation contained four vortices as opposed to the 
two displayed in Fig. |l[ This is merely a point of technical convenience. 
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Having put the gauge configurations on the lattice, the question can be ad- 
dressed what effect the center gauge fixing procedure has on these configu- 
rations. For the case of two thin vortices, fixing the configuration depicted 
in Fig. |I] (a) to maximal center gauge can be carried out trivially without 
any recourse to numerical computation. All one needs are the following two 
observations: The configuration shown in Fig. [1] (b) is gauge equivalent to 
the one shown in Fig. [I] (a) since all Wilson loops take the same values in 
the two configurations. At the same time, in Fig. [I] (b) all links take values 
in the center of the gauge group. Therefore, this configuration realizes the 
maximal possible value of the gauge fixing functional (0). Thus, the config- 
uration shown in Fig. [I] (b) is precisely a center gauge fixing image of the 
configuration shown in Fig. [l] (a). It is described by a string of links tak- 
ing a value corresponding to a nontrivial center element extending from the 
position of one thin vortex to the position of the other. This type of configu- 
ration, i.e. strings of links associated with a nontrivial center element of the 
gauge group, will in the following be termed a configuration of ideal vortices. 
It is, as just seen, gauge equivalent to a configuration of thin vortices at the 
positions where the strings end. 

It should be noted that there exist many other degenerate maxima of the 
gauge fixing functional on the same gauge orbit, namely all deformations of 
the string of links, cf. Fig. [I] (c). These are all images of the configura- 
tion shown in Fig. [I] (b) under gauge transformations which are elements of 
the center of the gauge group. This is not surprising; by construction, the 
gauge fixing condition (|6]) is invariant under such transformations (e.g. in the 
SU(2) case, it is insensitive to the sign of the link variables). In other words, 
it leaves the center of the gauge group unfixed. As a consequence, the precise 
trajectory of the string of nontrivial links describing ideal vortices is arbi- 
trary; only its two endpoints (corresponding here to the positions of the two 
thin vortices) are fixed, and they embody all gauge-invariant information. 

Note that in three-dimensional space-time, the sets of nontrivial links describ- 
ing ideal vortices generically form two-dimensional sheets, whereas in four- 
dimensional space-time, they form three-dimensional volumes. The bound- 
aries of the sheets or volumes, respectively, describe the locations of the thin 
vortices to which the ideal vortex configuration is gauge equivalent. 

As a further aside, note that center projection leaves the center gauge fixed 
configuration (Fig. [I] (b)) unchanged and thus does not truncate any phys- 
ical information contained in it. By contrast, center projecting the original 
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configuration (Fig. [I] (a)) before gauge fixing completely removes the vortices 
and yields a trivial lattice with all links at unity. 

Continuing to the case of a thick vortex, it is clear that this case cannot be 
solved in a similarly trivial fashion. Here, the latticized gauge configurations 
can be brought into the maximal center gauge numerically by maximizing 
the gauge functional (|6|) under the action of arbitrary SU (2) gauge transfor- 
mations on the sites of the lattice. As described above, lattice tori containing 
four vortices were used. Two different cases were examined: Two thick vor- 
tices of one orientation together with two thin vortices of the other orientation 
(cf. Fig. Q (a)), and one thick vortex together with three thin vortices, again 
such that the total flux piercing the plane vanishes (cf. Fig. |2| (b)). In order 
to make the new effects arising in the case of a thick vortex more visible, a 
Fermi distribution 

f( r ) — 1 I M 2^1 

A ' 1 + exp(5(r - 4/)/Z) 1 ' 

where I denotes the lattice spacing, was chosen for the profile function f(r) 
(cf. eq. (|10D). This not only spreads out the flux of the vortex but actually 
localizes it on a ring of radius 41 concentric to the core of the vortex. 



ill 



■ 



■ 



(a) 



(b) 



Figure 2: Configurations containing thin vortices and thick vortices, the flux 
of the latter being concentrated on the displayed rings. These configurations 
can be transformed into the maximal center gauge numerically (see text). 
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The gauge equivalent configurations which maximize the center gauge fixing 
functional (§) will contain a certain set of links which lie in the "south- 
ern hemisphere" of the £77(2) group, i.e. nearer to the center element (— 1), 
whereas the other links remain closer to the center element (+1). Center pro- 
jection then associates a collection of ideal vortices, or, gauge-equivalently, 
thin vortices, with the original configuration. Carrying out this procedure 
numerically for the configurations depicted in Fig. ^|, one obtains thin vor- 
tices which approximate the locations of the original thick vortices; however, 
the precise location of the thin vortex within the profile of the thick vor- 
tex turns out to depend on the precise realization of the center gauge via 
the arbitrary function g in the gauge fixing functional (|6]). Specifically, for 
the configuration containing one thick vortex, Fig. ^| (b), it turns out that 
for the gauge fixing function g(t) = t 2 , the thin vortex is located concen- 
trically with the thick vortex, whereas for g(t) = 1 — 6(2 — t)(l — t 2 /4:) 2 , 
the thin vortex is located one lattice spacing away from the center of the 
thick vortex. This latter effect becomes more pronounced if one uses the 
configuration with two thick vortices, Fig. ^ (a), again with the gauge fixing 
function g(t) = 1 — 6(2 — t)(l— 1 2 /4) 2 . There, the thin vortex is located four 
lattice spacings away from the center of the thick vortex, i.e. on the ring of 
flux depicted in Fig. ^| (a). 

The conclusion reached from these examples is therefore that the maximal 
center gauge does yield thin vortices which roughly approximate the locations 
of thick vortices in gauge field configurations. This corroborates the findings 



of pi , p5| ; note however recent indications of considerable instability in the 



numerical gauge-fixing procedure [19]. Furthermore, the precise coordinates 



of the thin center projection vortices turn out to be gauge-dependent, via the 
gauge fixing function g. Also, comparing the results obtained with the two 
different configurations depicted in Fig. ^ the precise location of a thin vortex 
approximating a given thick vortex evidently also depends on the structure of 
other vortices in the vicinity. Therefore, it seems that the center gauge fixing 
procedure introduces nonlocal interactions into the center projection vortex 
effective theory independent of any reference to the Yang-Mills dynamics 
which ultimately governs the original thick vortex configurations. 



2.3 Separating center and coset degrees of freedom 

In the attempt to formulate an analogue of the center gauge and center pro- 
jection procedure in continuum Yang-Mills theory, one obviously encounters 
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a conceptual obstacle: Continuum Yang-Mills theory is described in terms 
of a gauge potential A a ^(x) which is defined in the adjoint representation of 
the gauge group. Thus, in the continuum, the gauge symmetry of Yang-Mills 
theory is only a local SU(N)/Z(N), i.e. a coset, symmetry. Related con- 
siderations concerning this issue in the Schrodinger picture can be found in 



[2?j] . The first step in dealing with this problem is to carefully separate center 
and coset degrees of freedom in the lattice gauge fixing procedure discussed 
above. The treatment below will for simplicity mainly concentrate on the 
case of an SU{2) gauge group; however, the generalization to any SU(N) 
group is conceptually, if not necessarily notationally, straightforward. One 
only has to keep in mind that higher SU(N) groups have more center ele- 
ments, and that therefore there are more species of vortices, corresponding 
to the different possible fluxes. The only qualitative difference in terms of 
the topology of vortex configurations is that, for higher gauge groups (and 
three or more space-time dimensions), vortices may branch into several vor- 
tices with different center fluxes, as long as flux conservation is respected at 
branchings. 

To have a slightly more precise terminology at hand in the discussion below, 
it is useful to decompose the link variables U^(x) describing a SU{2) lattice 
configuration into a center part 

Z^(x) = sign tr U^(x) (13) 

and a coset part U' (x) defined by 

U,(x) = Z,(x)U'^x) . (14) 

The lattice described by Z^ix) will be referred to as the center 'part of the 
lattice configuration U^x), whereas the lattice described by U'(x) will be 
referred to as the coset part. Note that the center gauge fixing functional (^) 
depends only on the coset part of a configuration; center projection simply 
means mapping U^(x) onto Z^ix). 

Furthermore, below the terminology coset transformation will be used to 
characterize gauge transformations, cf. eq. ([!]), which change only the coset 
part of a lattice configuration. Note that this does not imply that the cor- 
responding unitary transformation matrices acting at lattice sites are neces- 
sarily from the northern hemisphere of the SU(2) group. Note furthermore 
that the set of possible coset transformations depends on the initial gauge 
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configuration and that in general not all members of this set can be contin- 
uously deformed to the identity while remaining in the set. Also, smooth 
gauge transformations G in the continuum formulation of Yang- Mills theory, 

A,^A° = G(A^ + d^ (15) 

are by construction coset gauge transformations. 

Clearly, if one puts a smooth continuum gauge configuration A^x) onto 
a sufficiently fine lattice, the associated center part is trivial, Z^{x) = 1. 
However, even in such a case, gauge fixing in general, and the center gauge 
fixing procedure in particular, may induce a nontrivial center part in the 
gauge fixed lattice; explicit examples were given in the previous section. In 
other words, some links in the gauge fixed lattice end up in the southern 
hemisphere, i.e. closer to the nontrivial center element (—1) than to the 
trivial element (+1); ideal vortices are introduced. 

In the 577(2) lattice description of the center gauge fixing procedure, one 
allows arbitrary SU(2) gauge transformations when maximizing the gauge 
fixing functional (|6]). This procedure can be broken down into smaller steps 
as follows. 

1. Start from a configuration with a trivial center part, i.e. with all links 
in the northern hemisphere, such as one obtains by putting any smooth 
continuum gauge configuration onto a sufficiently fine lattice. Apply a 
gauge transformation which rotates a certain arbitrary, but fixed, set 
of links to the southern hemisphere, i.e. such that the center part of 
the transformed lattice displays the nontrivial value (—1) on that set 
of linksQ 

2. Maximize the gauge fixing functional using only coset transformations, 
i.e. keeping the center part of the configuration invariant. In other 
words, links are not allowed to switch "hemispheres" of the SU(2) 
group. Note that this second step operates purely on the coset part 
of the lattice configuration obtained after step 1.) above; not only are 

7 As long as the original links are sufficiently close to unity (which one can always achieve 
by choosing the lattice spacing accordingly), this transformation can e.g. be effected by 
applying the group element V = exp(27rT3/3) to the site from which the links in question 
emanate and the group element to the site at which they end, i.e. U — > VUV, where 
U denotes the original link variable associated with the link under consideration. One can 
also exchange V and V', which may become necessary if two links in question share a site. 
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only coset transformations allowed, but also the gauge fixing functional 
@ depends exclusively on the coset part of the link variables. 

3. Finally, this two step procedure must be repeated for all possible first 
step choices of the center part of the transformed lattice, in order to 
cover all possible gauge transformations contained in the full SU(2) 
gauge group. For each of these choices, corresponding to a certain 
ideal vortex configuration, one obtains in the associated second step a 
constrained maximum of the gauge fixing functional @; from all these 
maxima, one must choose the highest one in order to obtain the gauge- 
fixed configuration reached under the action of the full SU (2) gauge 
group. 

Slightly rephrased, the procedure described above splits the process of cen- 
ter gauge fixing under the SU{2) gauge group up into the following parts: 
One introduces, by a suitable gauge transformation, into an initially smooth 
gauge configuration ideal vortices at fixed, but arbitrary, locations. This step 
transfers information from the coset part into the center part of the config- 
uration. From this point onwards, however, center and coset parts remain 
strictly separated and subsequent gauge transformations operate only on the 
coset part. The latter is transformed, by using purely coset transformations, 
to the (constrained) maximum of the (purely coset-dependent) gauge fixing 
functional (||). The algorithm must be repeated for all possible choices of 
ideal vortex configurations. The choice which leads to the global maximum 
of the gauge fixing functional determines the center gauge and center projec- 
tion image of the original configuration. 

Formally, one can summarize this as follows. The center gauge fixing condi- 
tion, cf. eq. (|B|), corresponds to the two-fold maximization problem 

(16) 

where U^x) is the gauge configuration one starts with and Z denotes the 
center part of the transformed lattice chosen in step 1.) above, realized by 
the gauge transformation V[Z] as specified there. The prime emphasizes ex- 
plicitly that only the coset part of the resulting configuration is considered 
further, although this is already guaranteed by the absolute value prescrip- 
tion. In other words, the configurations {U v ) f are stripped of all ideal vor- 
tices. Finally, G denotes all coset gauge transformations considered in step 
2.) above. 



max max q 
z g ^ 



lG 
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3 Taking the continuum limit 



3.1 Continuum limit of lattice configurations 

As the first step in formulating a continuum version of the maximal center 
gauge, it is necessary to define the continuum analogues of the various lattice 
configurations considered above. 

3.1.1 Continuum limit of the coset part 

The coset part of a lattice configuration can be directly associated with a 
regular continuum gauge field constructed such as to reproduce all links of 
the coset configuration if put on the lattice. One could e.g. achieve this con- 
struction by giving an "inverse blocking" prescription, successively replacing 
coarser by finer lattices while dividing up the phases carried by the coarse 
links onto the fine links such as to achieve a smooth continuum limit. In this 
limit, it will be possible to expand an exponential representation of the link 
variables 

Ufa) = exp{-Ufa)) (17) 

into powers of the lattice spacing / and regard A'Jx) as the continuum gauge 
field. As a side remark, note that if one wished to apply such an "inverse 
blocking" procedure to the center part of a configuration, the corresponding 
prescription would have to be radically different in order to preserve the 
center, or ideal vortex, character of the configuration: If a nontrivial center 
link of a coarse lattice is divided up into two links of a finer lattice, the 
whole nontrivial center phase would have to be associated with one of the 
two links, whereas the other link would be associated with unity. The phase 
is not divided up among the two links, in contrast to the coset case, and 
the continuum limit is singular. In particular, expanding an exponential 
representation analogous to ( |TTD is by construction never possible for the 
center part of a configuration. Due to these difficulties, a different approach 
will have to be taken when discussing the continuum limit of the center part 
below. 

Before proceeding, it is worth noting that the gauge freedom associated with 
the continuum gauge field A'^ in (|T7| ) matches the one associated with the 
coset lattice U' in step 2.) of the center gauge fixing procedure as described 
in section |2.3| ; namely, it consists of all possible coset transformations. Thus, 
in describing a coset lattice in the continuum limit by means of a gauge 
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field A'^ as in flT7|), one is not introducing any restriction on the residual 
gauge freedom allowed in step 2.) of the center gauge fixing procedure. The 
additional gauge freedom implied by the full SU(N) symmetry present in 
lattice Yang-Mills theory has been explicitly separated off into the choice of 
ideal vortex configuration in step 1.). 

3.1.2 Continuum limit of the center part 

The center part of a configuration contains exclusively ideal vortices, de- 
scribed by sets of links taking a nontrivial value from the center of the gauge 
group. At this point, it is useful to give an alternative characterization of 
such a configuration which does not make explicit reference to a space-time 
lattice. This can be achieved via the effect the configuration has on arbitrary 
Wilson lines. Consider therefore calculating an arbitrary Wilson line; as one 
travels along the line, gathering up phases from the different links making 
up the path, the Wilson line either stays constant or picks up a phase cor- 
responding to a nontrivial element from the center of the gauge group on a 
length of one lattice spacing /. Expressed through the vortex content, this 
latter change happens whenever the Wilson line intersects a string, surface, 
or volume of nontrivial links (depending on the space-time dimension) de- 
scribing an ideal vortex configuration. As the lattice spacing is taken to 
zero, these strings, surfaces, or volumes become infinitely thin. Therefore, 
the continuum analogue of an ideal vortex configuration consists in specify- 
ing strings, surfaces, or volumes E (depending on the space-time dimension) 
which, when intersected by a Wilson line, contribute a center element factor 
to the latter. 

3.2 Singular gauge transformations 

In the previous sections, the effect of the maximal center gauge on vortex 
gauge fields originally defined in the continuum theory was investigated by 
first putting the Yang-Mills field on the lattice, subsequently implementing 
the maximal center gauge condition (which, so far, is defined only on the 
lattice) and afterwards considering again the continuum limit. One would, 
however, prefer not to use the detour via the lattice, but rather perform 
the maximal center gauge fixing directly in the continuum. Particular care 
is needed in defining the analogue of step 1.) of the maximal center gauge 
fixing procedure as laid out in section |2.3| . It should be clear from the previous 
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discussion that this step in the continuum will introduce singularities (the 
ideal vortices) into the initially smooth gauge configurations. As already 
partly indicated above, these singularities are to be defined via their lattice 
origin, taking care in particular that gauge invariant quantities are unaffected 
by gauge transformations, even when the latter have a singular continuum 
limit. Keeping this methodology in mind, it will be possible to paraphrase 
the center gauge fixing procedure, as described step by step in section [2~3l , 
directly in the continuum. 

Consider a continuum Abelian gauge transformation V(k, S, x) specified as 
follows. Let V be continuous everywhere except at hypersurfaces S of di- 
mension D — 1 if D is the dimension of space-time. At such hypersurfaces, V 
shall have the following property. Let t denote the local coordinate perpen- 
dicular to E, with t = on S, and let x± denote the other local space-time 
coordinates. Then 



V(E,x ± ,t=-e)d^(E,x ± ,t=-e) = V(E,x x ,t = e)d i y\E,x ± ,t = e) 



where e — > and k labels the N — 1 different nontrivial center elements Z(k) 
of the SU(N) gauge group under consideration. In other words, V jumps 
by a center phase Z(k) at hypersurfaces S. An explicit realization of such a 
gauge transformation V will be given in the next section. 

Consider applying a gauge transformation such as characterized above to an 
initial smooth field configuration A, 



The first term on the right hand side still represents a smooth field everywhere 
except possibly on the edge <9£ of S, since it is quadratic in V and hence 
insensitive to the center phase V picks up at hypersurfaces S. More care 
must be taken in properly defining the singularity in the second term. For 
this purpose, consider a Wilson line describing a straight-line path C from 
(x±, t = — e) to (x±, t = e), with x± and t defined as above, and infinitesimal 
e. Since the Wilson line transforms as 



W[A]{C) -> W[A V ]{C) = V(x ±1 t= e)W[A}{C)V\x^t = -e) (20) 



V(k, S, xl, t = e)V\k, S, x ±1 t = -e) 



Z(k) 



(18) 



A -> A v = VAV^ + VdV ] 



(19) 



and furthermore W[A](C) = 1 + 0(e), one has 

W[A v ](C) = Z(k) + 0(e) 



(21) 
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Therefore, the singularities of A v at the hypersurfaces E describe nothing 
but an ideal vortex configuration A(T,). On the other hand, the term VdV' 
is a pure gauge. Therefore, the regular part of VdV^ must compensate the 
ideal vortex part described by the singularities at E such as to render all 
(closed) Wilson loops taken in the configuration VdV^ trivial. This means 
that VdV^ can be decomposed asQ 

V{T)dV\T) = A(Z) -a(9S) (22) 

with the ideal vortex configuration A(T0) localized on a hypersurface S, and 
a describing thin vortices localized on <9E. The thin vortex field a can be 
defined on E as 

a|s = — lim V(E, x±, t = — e)dV j! ('E, x±, t = — e) 

e— >0 

= - lim V{E, x ± , t = e)dV*(E, x ± , t = e) , (23) 
cf. eqs. ( fL8|) . Away from E, 

a = -VdV ] (24) 

since the ideal vortex configuration A(H) only has support on E. 

Before applying these singular gauge transformations to implement the max- 
imal center gauge in continuum Yang-Mills theory, a few further properties 
should be noted for convenience. Given a general thin vortex configuration 
a, a corresponding gauge transformation V generating it can always be re- 
covered via 

V(y) = exp (j^ ^ a^dx^ (25) 

using an arbitrary fixed starting point y^. Path ordering is unnecessary, since 
a was constructed to be an Abelian field. Besides the trivial ambiguity in the 
choice of yo, which merely leaves a constant phase undetermined in V, there 
is a more important ambiguity, namely how the path P(yo, y) circumvents the 
locations of the thin vortices contained in a. This choice precisely determines 
the hypersurfaces E at which V jumps and which therefore describe the 
associated ideal vortex configuration A. 

From this, it also becomes clear that the thin vortex configuration a can 
be chosen independently of the precise singular surfaces E describing the 

8 The minus sign in front of a(9S) is a convention introduced for later notational 
transparency. 
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associated ideal vortex configuration. Merely the boundary <9E must coin- 
cide with the location of the thin vortices contained in a. Indeed, given a 
transformation Vi(Ei) inducing a thin vortex field a, one can explicitly give 
another transformation I4(E 2 ), with <9£i = <9E 2 , generating the same field a 
as follows. Let M denote the interior of the closed hypersurface[] Si U — E 2 . 
In D space-time dimensions, M is D-dimensional. Then define 

This transformation generates the same thin vortex field a as Vi(Si) and an 
ideal vortex configuration described by the hypersurface E 2 . In particular, 
using Vi(Si) and V 2 (E 2 ), one can deform the ideal vortex part of an arbitrary 
gauge configuration without changing the configuration in any other respect: 

(A + A^))^ 2 =A + A{E 2 ) (27) 

This is the continuum analogue of a transformation from the center of the 
gauge group. It exclusively deforms ideal vortices. Furthermore, the form of 
( |2"o| ) shows that a homogeneous rotation of any gauge field is independent of 
the location of the singular hypersurfaces S, 

K(E 1 )Ay 1 t (E 1 ) = ^ 2 (E 2 )Ay 2 t (E 2 ) (28) 

as long as V\ and generate the same thin vortex field a, i.e. in particular 
dEi = <9E 2 . 

As a simple example, consider the SU(2) transformation 

V = exp(y?T 3 ) (29) 

in the plane described by the polar coordinates r, (p with (p e [0, 27r[. At 
<f — 0, one has 

V(r, if = e)V f (r, tp = 2tt - e) = -1 + 0(e) (30) 
and, away from tp = 0, 



a r 



--Vd v V* = ^T 3 (31) 
-Vd r V ] = (32) 

9 The notation "— S" is introduced to specify the orientation of a hypersurface. The 
hypersurface — S is located on the same set of space-time points as the hypersurface S, 
but is associated with the reverse orientation. 
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i.e. one observes an ideal vortex at cp — and a thin vortex (cf. eq. ([])) at 
the origin, r = 0. By simply changing the interval on which (p is defined (in 
general, as a function of r), one can arbitrarily deform the ideal vortex A(T,) 
without changing a(<9£). 



3.3 Continuum representation of center vortices 

In the continuum limit, ideal vortex configurations Ap(k, S, x) are specified 
by hypersurfaces £ (strings, sheets, or volumes in D = 2, 3 and 4, respec- 
tively) which, when intersected by a Wilson line C, contribute a center ele- 
ment Z(k) to the latter. 

An explicit (singular) gauge field representation of an ideal vortex configu- 
ration in D space-time dimensions is given by 

A^x, T)=E I d D ~ 1 a ll 8 D {x - x{a)) . (33) 

Here, £ describes the infinitely thin continuum limit of the string, surface, or 
volume of nontrivial center element links making up the ideal vortex config- 
uration. £ is D — 1-dimensional and its boundary S = <9£ gives the location 
of the thin vortex a^(<9£, x) to which the ideal vortex configuration *4. M (£, x) 
is gauge equivalent (see below). Furthermore, x^{a) = x^ai, 02, . . . , &d-i) 
denotes a parametrization of the (D — l)-dimensional hypersurface S. Also, 

d CT^i — j-yj e^tof!. ..a D -i ^ "'ai-ao-i (34) 

is the dual of the (D — l)-dimensional volume element 

iD—l 9x ai dx aD1 foc\ 
d V ai ...a n -i — e k x ...k D -i-R ~r dO\ . . . da r>-\ ■ (35) 

Note that these volume elements also define an orientation of the vortex. 
Like in the special case (|D|), in general two different orientations of vortices 
can be distinguished. 

On the other hand, the color structure of the vortex is encoded in E. Note 
that the following discussion of the color factors E applies equally to thin 
vortices, to be treated in detail further below. The different kinds of vortex 
configurations introduced in this work are distinguished from one another 
only by their space-time structure, but they carry identical color structure. 
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It should furthermore be emphasized that the orientation of vortex flux is 
encoded already in the space-time structure, as mentioned above; this entails 
consequences regarding the set of color factors E necessary to generate all 
possible distinct vortex fluxes. For example, the substitution of E by — E 
should not be considered to lead to a new class of vortices; the configurations 
generated by this substitution are already accounted for by keeping the orig- 
inal color structure E, but reversing the orientation of the hypersurface S 
in (|3lf). As specified below, E should be restricted to a certain fundamental 
domain. 

For the gauge group SU(N), ideal center vortices in the continuum are given 
by eq. (|33|) with E(k) = E c (k)T c satisfying 

e ~ E ( k ) = Z{k) , k = l,...,N-l (36) 

where only generators T c from the Cartan subalgebra are associated with non- 
vanishing coefficients E c (k), and Z(k) denotes the N — 1 nontrivial center el- 
ements of the SU(N) group. The E c (k) can be considered as the components 
of vectors E(k) in the Cartan subalgebra which, up to a factor 2tt, coincide 
with the so-called co- weights i.e. E(k) = 2tt fi(k). The co-weights f/,(k) 
are dual to the simple roots and define the corners of the fundamental do- 
main of the SU(N) algebra. They are defined in the adjoint representation 
{T a )bc = fabc by 

e~ & = 1 , E = E c f c (37) 

Here, mainly the gauge group SU(2) will be considered, for which only one 
nontrivial center element Z — — 1 exists, associated with the color structure 

E = E 3 T 3 , E 3 = 2tt . (38) 

For the gauge group SU(3), there are two nontrivial center elements 

Z(l) = e i27r/3 , Z(2) = e- i27r/3 (39) 

which are generated from (|36| ) by the vectors E(l) and E(2), whose non- 
vanishing components explicitly read 

47T 2tt 
£ 8 (1) = _, £ 3 (2) = 2vr, E 8 (2) = -j= . (40) 

Note that both vectors E(l) and E(2) satisfy the normalization 

tr(E(k)) 2 = - l -E c {k)E c {k) = (41) 
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while 



4vr 2 

tr£(l)£(2) = -— . (42) 



With (pop , one has for the Wilson line along a path C in the vortex configu- 
ration (|33|) 

W[A]{C) = Z(k) I(c '^ (43) 

where 

/(C,E) = / ^ / d^cy^x-ac^)) (44) 

is the intersection number between C and E in D dimensions. For closed 
paths C, i.e. for Wilson loops, this equals the linking number between C and 
S = dE, 

I{C, E) = L{C, S = dE) . (45) 

In particular, A therefore satisfies the properties characterizing an ideal 
vortex configuration as formulated at the beginning of this section. As al- 



ready mentioned in connection with eq. (27), the sheet E can be arbitrarily 
deformed by a (singular) Abelian gauge transformation, while keeping its 
boundary S = <9E fixed; this is entirely analogous to the lattice, where the 
stack of links U — — 1 (in the case of SU(2)) forming the hypersurface E 
attached to the center vortex S = <9E can be deformed by center gauge 
transformations. A more explicit construction of this transformation is given 
further below in eq. fl65|). 

On the other hand, in order to find a general representation of thin vortex 
configurations, consider the following gauge transformation, 

V(k,E,x) = exp(E(k)n(E,x)) (46) 

where 

n(E,x) = —l dP~ x a % (47) 



oo D ^i Js ( x ( a ) — X ) D 

is the solid angle taken up by the hypersurface E when viewed from x. Here, 

2tt d / 2 

= f(W) (48) 

is the area of the unit sphere Sd-i in D dimensions, so that the solid angle as 
defined above is normalized to unity for a point x inside a closed surface E. 
Note that the solid angle is defined with a sign depending on the orientation 
of E as rays emanating from x pierce E. A deformation of E keeping its 
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boundary <9E fixed leaves f2(E, x) invariant unless x crosses E. In the latter 
case, f2(E,x) changes by an integer. 

For later convenience, note also that the solid angle fl47|) can be represented 

as 

fi(E,x) = f d D - l a^D(x - x{a)) (49) 

where D(x) is the Green's function of the D-dimensional Laplacian, 

- d*d*D(x - x') = 5\x - x') . (50) 

For a closed surface E = dM the solid angle agrees with the characteristic 
function x(M, x) of M by Gaufi' theorem, 

tt(E,x)= [ d D a5 D (x -x(a)) = X (M,x) (51) 

J M 

i.e. any point inside E = dM is associated with a solid angle 1, any point 
outside E with a vanishing solid angle, 

^ I ' = {jothe^ < 52 > 

In this respect, the normalized solid angle may also be interpreted as the 
linking number of the point x and the hypersurface E. 

Consider now the gauge field induced by the gauge transformation (f|6]) de- 
fined by the solid angle f2(E,x), 

V(k, E, x)d i y\k, E, x) = -Eifyd^nCE, x) . (53) 

Obviously, one has 

I dx^ d^{E, x) = (54) 

J C 

for any closed curve C, so that any Wilson loop receives the trivial contribu- 
tion 

W[VdV ] }(C) = 1 . (55) 
Now one can show that indeed 

V{k, E, x)d fM V\k, E, x) = A„(k, E, x) - a^k, <9E, x) (56) 

where a M and A^ are the thin and the ideal vortex fields, respectively, de- 
scribing the same magnetic flux localized on <9E. While A^ has already been 
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defined by eq. (|3"3"|), an explicit representation for the thin vortex a^k, <9£, x) 
in ( |56|) will be derived below. For this purpose, consider 

M fi(£, x) = j d z ~a v d^d v D(x - x(a)) (57) 

and write 

d„d v = 5^d 2 - (d^d 2 - dM (58) 

i.e. express the longitudinal projector L^ v = d^d u /d 2 as 5^ u — V^ u , where 
Vfiv = — d^dy/d 2 is the transversal projector. Inserting this relation into 
eq. (0) and using the definition of the Green's function fl5Dp, one obtains 

d li n(£,x)=-J ^^/(i-s^J-J d^cr^S^-d^Dix-xia)) . 

(59) 

Upon inserting this relation into (|53"D, the first term produces precisely the 
ideal vortex field A(k, S, x) as defined in eq. (|33]) . To show that the second 
term yields the thin vortex field, it is necessary to express the transversal 
differential operator as 

5^ u d 2 - d^d v = ^e^ Ka pe„xapd K d x . (60) 

After straightforward calculation, using Stokes' theorem, one eventually ob- 
tains 

-/ d D - x a v (8^&-dMD{x-x(<r)) = - [ d D ' 2 ~a, K dlD{x-x{a)) (61) 

JT, JdT, 

Inserting this into eq. (|53|) , one finds the thin vortex for arbitrary <9E, 

a^k, dZ, x) = -E(k) [ d D - 2 d^d s K D(x - x(a)) (62) 

In Appendix |A.2| , it is shown explicitly that a^k, dH,x) in eq. (|62| ) indeed 
represents a center vortex, namely 

I dx^a^k, <9£, x) = E(k)L(C, dE) (63) 
J c 

where L(C, <9£) is the linking number between C and <9E. 

Note that, unlike the ideal vortex A(E,x), the thin vortex a(dTi, x) depends 
only on the boundary <9£, where the magnetic flux associated with the vortex 
is located. 
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It is now easy to see that a gauge transformation which deforms the singular 
sheet Ei, describing an ideal vortex, into another sheet E 2 while keeping the 
boundary <9E fixed is given by 

V^Ei^x) = V- l {k^x)V{k^ 2 ,x) (64) 
= exp [-£(&) (Q(Ei,x) -ft(£ 2 , x))\ (65) 

Since the boundaries of Ei and E 2 match, dH\ = <9E 2 , the thin vortices 
arising from V A (Ei)(9V A t(Ei) and V(E 2 )<9V^(E 2 ) cancel and one obtains 

(A(k, Ei, x ))V{k,^,x) = A{k ^ E2) x) _ (66) 



One can also easily convince oneself that the transformation fl65|) agrees with 
the one introduced in section |372| eq. (fffi). Since EiU(— E 2 ) = dM represents 
the closed orientable surface of some volume M, application of GauB' theorem 
yields 

n(Zi,x)-n(Z 2 ,x) = - I dP^d^Dix - x(a)) (67) 

+ / d D - x ap&*D{x - x(a)) 
= - [ d D - l audlD(x - x(a)) (68) 

•/£iU(-£ 2 ) 

= / d D a{-d?X)D{x-x{a)) . (69) 



Using the definition of the Green's function of the D-dimensional Laplacian 
(|50D, one obtains 

n(E u x) -n(E 2 ,x) = [ d D a5 D (x-x(a)) = X (M,x) (70) 

where x{^i x ) is the characteristic function of M defined in eq. (|H), and 
dM — Ei U (— E 2 ). Hence, the gauge transformation ( |65D can be expressed 
as 

Vik^^x) = exp[-E(k) X (M,x)} (71) 
= Z(k) x{M > x) . (72) 

Thus, one recovers precisely the transformation already introduced in 
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3.4 Continuum limit of the Maximal Center gauge 

On the basis of the groundwork laid out in the previous sections, it is now 
straightforward to construct the analogue of the maximal center gauge for 
continuum Yang-Mills fields. The transformations V(E) characterized by 
fll8|) allow to carry out precisely step 1.) of the maximal center gauge fixing 
procedure, as specified in section |2.3| , in the continuum. Given an initial 
smooth gauge field A, the transformed field A v ^ contains ideal vortices 
described by hypersurfaces £ of one's choice. Furthermore, in analogy to the 
lattice, where taking the coset part of a configuration amounts to leaving 
away all ideal vortices A(k, £, x), one can define the stripped configuration 



explicit representation being given by eq. fl6"2"|). This stripped configuration 
induces the same behavior for arbitrary Wilson line integrals as the config- 
uration A v , except for the center element jumps at the hypersurfaces E. In 
fact, (A v )' depends only on the thin vortex configuration a, with thin vor- 
tices at <9£, and not on the specific choice of the entire hypersurface E, cf. 
eq. (^8|) . It is therefore the precise analogue of the coset part in the lattice 
case. Note that {A v )' is still singular on <9£. 

According to step 2.) of the maximal center gauge fixing procedure as speci- 
fied in section [2.3| , one must now consider all coset gauge transformations of 
{A v )' such as to maximize the gauge fixing functional. As explained in the 
previous section, this step can be straightforwardly translated to the contin- 
uum, where gauge fields are subject to arbitrary local coset transformations. 
Expressing therefore the coset link configuration in the continuum limit, cf. 




(73) 



with the thin vortex field a defined 




eq. (0), 




(74) 



the maximal center gauge fixing condition flTBj ) can be rewritten as 




(75) 



where the arbitrary gauge fixing function has also been redefined, 




(76) 
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Like g, g should be a monotonously rising function. Note that in (173), the 
maximization is, as written, over all hypersurfaces E, and not over all possible 
associated transformations V(Tj). For every E, there are many possible V(H). 
Only one ^(E) must be constructed and inserted in ([Tj|) for every choice of 
E. 

Inserting (|73|) , the gauge fixing condition can be written as 

max max J d D x g (tr (GV{A - a)V^G ] + GdG^ . (77) 

This expression can be further simplified such as to eliminate all explicit ref- 
erence to the singular transformations V. In the following calculation, there 
is a slight subtlety. Namely, the integrand in ( |77| ) is continuous at hyper- 
surfaces E; it was precisely constructed to be so. Thus, one may replace its 
value on hypersurfaces E by the limit as one approaches E. In the following, 
the integrand in ([F_7]) will therefore be evaluated outside of E and its value 
on E is then defined by continuity. In particular, this means that one may 
set VdV* = -a, cf. eq. (||). 

By writing G = VGV^ , the argument of the gauge function g in (|77|) becomes 

tr (VG{A - a)G ] V ] + VG(V ] dV)G ] V ] + V(Gd&)V* + VdV^f (78) 

= tr (VGAGW + V(Gd&)V ] + VdV^f (79) 

= tr (VA°V^ -of (80) 

= ti(A Q -af (81) 

where the fact that V and a commute was used. Note that G is continuous 
on E, since it is quadratic in V and hence insensitive to the center phase V 
picks up at hypersurfaces E. Thus, it can still be considered a coset gauge 
transformation. Vice versa, any coset transformation G is associated with a 
coset transformation G. Therefore, the gauge fixing condition ([771) can be 
rewritten as 

max max J d D x g (tr [A a - a)^ (82) 

with the coset maximization running over all coset transformations G. 

In this form, maximal center gauge fixing is revealed more clearly as an 
approximation problem in the mathematical sense; it corresponds to approx- 
imating the gauge orbit of a given gauge configuration A as well as possible 
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by thin vortex configurations a. The precise metric defining how close a thin 
vortex field a is to the gauge orbit of A is encoded in the arbitrary function 
g. Thus, the continuum formulation of the maximal center gauge provides 
a deeper understanding of the empirical fact drawn from lattice experiments 
[fL0|l , namely that center projection vortices provide a rough localization of 
thick vortex structures present in full lattice Yang-Mills field configurations. 
Note that, in mathematical terms, the space of thin vortex fields a has no 
scalar product or orthogonality properties. Eq. (^) specifies a norm which 
serves to define an approximation problem within this space, namely maximal 
center gauge fixing. 

Note that the gauge fixing functional (|8^) only depends on the thin vortex 
configuration a, with vortices localized on <9£, and not on the whole choice 
of S. This is the continuum counterpart of the lattice observation that the 
maximal center gauge leaves the center of the gauge group unfixed, i.e. does 
not depend on the center part of a configuration. 

As a last step, it is interesting to recast the coset gauge fixing condition in 
differential form as follows. Consider £ and an associated thin vortex field 
a(<9£) as given. Then a gauge field A satisfies the condition of providing the 
maximal value of the gauge fixing functional under the action of arbitrary 
coset transformations G if the value of the gauge fixing functional (^) is 
stationary under infinitesimal transformations G = exp(—9) « 1 — 9, 

= — — [ d D xg(ti([A,9] +09 -a) 2 ) (83) 



69*(y) 

J d D x g' (tx{A - a) 2 ) {A c - a c ){-85 ac - A d f dac )5 D (x - y) (84) 
' (tr(A - a) 2 ) {Al - o»)] - g' (tr(A - a) 2 ) f adc A d a c (85) 
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where the prime on g denotes the derivative with respect to the argument. 
By multiplying with generators T a , this reduces to 

= F[A,a] = 2g" (ti(A-af) ■ ti ({A v - a v )d^A u - a u )) ■ (A^ - a ^) 

+g' (tr (A - a) 2 ) + a M , A„] - d^) (86) 

In the simple case g(t) = t, the gauge condition simplifies to 

[d„ + a M , A,j\ - <9 M a M = (87) 

If one in particular chooses the thin vortex configuration to satisfy the Landau 
gauge, d^a^ = 0, then this condition is nothing but the background gauge 
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with the background given by the thin vortex field a. Note that the thin 
vortex configuration fl62|) satisfies the Landau gauge, as can be inferred from 
eq. ©. ' 

Furthermore, if one has picked out of the gauge orbit of a given configuration 



A the element A [a] which satisfies (|86|) for fixed arbitrary a (coset gauge 



fixing), then the remaining part of the maximal center gauge fixing procedure, 
i.e. finding the optimal a, reduces to solving 

max J d D x g (tx(A[a] - a) 2 ) . (88) 



As mentioned repeatedly above, all choices of E in (pq) corresponding to the 
same <9E are degenerate due to the unfixed center part of the SU{2) gauge 
freedom. Also, for every <9E, it is only necessary to construct one particular 
a(<9E), cf. the discussion after eq. (|76|) . 

Finally, having obtained an optimal hypersurface E, center projection simply 
means replacing the full gauge field by the corresponding ideal vortex con- 
figuration -4(E), cf. eq. (^). This center projected configuration can then 
e.g. be used to evaluate observables such as the Wilson loop, cf. eq. P3]). 



3.5 Remarks on the center gauges in continuum Yang- 
Mills theory 

3.5.1 Revisiting the thick vortex 

In general, the gauge fixing condition (^) is not solvable by analytical means; 
at this stage, the value of the continuum formulation of the maximal center 
gauge lies more on a conceptual than on a practical level. However, it is 
worthwhile to briefly reexamine some properties of the examples presented 
in section |2.2| in the lattice framework. Namely, consider approximating (the 
gauge orbit of) a given thick SU(2) vortex configuration centered on the 



origin of a two-dimensional space-time plane, cf. eq. ( 10) 



A = T, 



f{^WTy I ) 



x 2 + y 2 

by a thin vortex centered at the point xoe x 

1 



(ye x - xe y ) , dA = 



a = Tk 



[x - x ) 2 + y'- 



[x — xo)e 



yj 



(89) 



(90) 



30 



where the offset xq is still to be varied; e^, e y denote the Cartesian unit 
vectors in the plane. For definiteness, the following treatment was carried 
out using a profile function f(r) = 6(r — R), i.e. a ring-shaped thick vortex 
of radius R similar to the one considered in section |2.2| . In the case Xq = 0, 
the gauge condition (|36|) is satisfied by these configurations, for arbitrary- 
gauge fixing function g. For xq 7^ 0, the condition fl86|) is in general not 
satisfied. Nevertheless, in terms of the variational form ( jS2|) of the gauge 
fixing condition, one has the relation 

J d D x g (tr(A - a) 2 ) < max J d D x g (ti(A° - a) 2 ) (91) 

where equality is guaranteed for x = 0. Now, it is possible to find gauge 
fixing functions g, e.g. g = — tanh(i? 4 t 2 ), for which the quantity 

J d D xg{ii(A-a) 2 ) (92) 

appearing on the left hand side of (0), viewed as a function of x , exhibits 
a minimum at xq = 0. In view of (|9l|) , this implies that also the quantity on 
the right hand side, 

max f d D x g (ii{A G - a) 2 ) (93) 

exhibits a minimum at Xq = (remember that equality is guaranteed in (|9T|) 
for x = 0). Therefore, according to the variational form of the gauge fixing 
condition (p2|), the optimal location of the approximating thin vortex con- 
figuration a, parametrized by xq, is achieved for xq 7^ 0. This confirms the 
phenomenon observed in the lattice examples in section |2.2j , that the thin 
vortex induced by applying the center gauge fixing procedure to a radially 
symmetric thick vortex configuration does not necessarily appear concentri- 
cally with the original thick vortex. This property thus does not constitute 
a lattice artefact, but persists also in the continuum theory. 

As a final remark, note that the simple choice of gauge fixing function g(t) = t 
leads to a trivial optimal thin vortex configuration a = for any smooth 
initial configuration A. This is due to the fact that a thin vortex field diverges 
as 1/r, where r is the distance from the vortex, and therefore tr(A — a) 2 
diverges as —1/r 2 at the positions of the thin vortices. Therefore, in the case 
g(t) = t, the only choice of thin vortex field which does not render the gauge 
fixing functional (|85p divergent is a = 0. Consequently, a = trivially turns 
out to be the optimal thin vortex approximation, for any smooth A^x). The 
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choice g(t) = t for the gauge fixing function thus at first sight does not seem 
very useful, since it cannot lead to a faithful representation of the vortex 
content of the original smooth gauge field A. The same is true for any g(t) 
which diverges as fast as, or faster than, g(t) — t as t — > oo. 



3.5.2 Remarks on the Direct Maximal Center Gauge 

The choice of gauge fixing function g(t) = t is the continuum analogue of a 



variety of lattice gauge fixing functions g(x), cf. eq. ( |76| ) in the limit / — > 0. 
Among these is the case g(x) = x 2 (after leaving away irrelevant constants). 
This is the standard "Direct Maximal Center Gauge" P1,P1,[^1- Strictly 



speaking, therefore, this gauge does not have a useful continuum limit on the 
level of gauge fixing individual smooth continuum gauge field configurations. 
The question arises how this gauge leads to a finite center projection vortex 
density in the continuum limit of the lattice formulation, as evidenced by the 
renormalization group scaling behavior of this density |20] (note erratum in 



1 1), see also [21|. Several possible explanations come to mind. 



For one, in a lattice calculation, one never takes the continuum limit of the 
field configurations at an intermediate stage. Instead, one operates with 
regularized lattice configurations and only extrapolates observables at the 
end of the calculation to the continuum limit. Thus, on the lattice, one 
always stays away a certain minimal distance from the 1/r-singularity of 
the continuum thin vortex, i.e. one never probes the t — > oo behavior of 
g(t). More precisely, one cannot make a distinction between g(t) = t and a 
"capped" version grit) = t9(T — t)+T9(t — T), with a sufficiently large value 
of T. Of course, as the lattice spacing is taken to zero, T will have to diverge if 
g and gx are to remain equivalent. Now, it is entirely possible that the choice 
g(t) = t leads to scaling violations of the projection vortex density on very 
fine lattices, such that this density ultimately tends to zero. On the other 
hand, a capped version grit) as above, with large but fixed T, allows a finite 
projection vortex density for smooth continuum configurations. Possibly, 
present lattices are still too coarse to reveal the above scaling violations and 
distinguish between g(t) = t and grit) with large fixed T. An observation 
which supports this possibility is that on present lattices, one still finds in 
typical configurations an abundance of negative plaquettes, i.e. structures 
indistinguishable from thin vortices on the scale of the lattice spacing. On the 
other hand, the appeal of this explanation is lessened by the fact that it calls 
into question the relevance of the scaling behavior observed for the projection 
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vortex density; usually, such behavior is considered a firm indication that the 
continuum limit of an observable can be extrapolated with confidence. 

A second explanation can be based on the fact that functional integrals in 
field theory are dominated, for entropy reasons, by configurations with infi- 
nite action, not by smooth configurations For instance, if typical con- 



figurations A contain thin vortices, despite the fact that they are associated 
with a divergent Yang-Mills action, then in the process of maximal center 
gauge fixing, i.e. approximating the configurations by thin vortices, one 
would indeed introduce a nontrivial thin vortex configuration a to cancel 
all thin vortices in the combination A — a. This would happen also for the 
choice of gauge fixing function g(t) = t. However, it has been argued that 
thin vortices are too singular to survive the continuum limit of lattice Yang- 
Mills theory f[? |. In this simple form, this second explanation is therefore 



questionable. 

However, a third explanation which is a slightly more sophisticated variant of 
the second one is possible. When fixing to the maximal center gauge on the 
lattice, one never finds the exact maximum of the gauge fixing functional, 
but only an approximation of this maximum. Put another way, consider 
introducing the (negative) gauge fixing functional as a weight into the Yang- 
Mills partition function, 

Sqf = ~qJ d D x g (tr(A[a] - a) 2 ) (94) 

cf. eq. (|88|). Then, in practice, one does not insist on the limit q —>■ oo, which 
would imply finding the exact minimum of Sgf', instead one may be satisfied 
with a finite value of q. In this competition between the gauge fixing 

term and the Yang-Mills action (and also a competition between different 
Gribov copies, i.e. relative minima of Sgf) becomes possible. Even in a 
case such as g(t) = t, where the gauge fixing term in general will diverge for 
nontrivial thin vortex configurations a, these divergences may be swamped 
by the infinite action of typical Yang-Mills field configurations. Note that for 
this to work, only Yang-Mills configurations considerably less singular than 
the thin vortices invoked in the second explanation above are necessary. E.g., 
in the case of the gauge fixing function g(t) = t, the gauge fixing term for 
a nontrivial thin vortex configuration a will diverge as / drr(l/r 2 ), where 
r denotes the distance from the thin vortex. This is independent of the 
dimension of space-time; the coordinates perpendicular to a vortex always 
span two dimensions. On the other hand, in four space-time dimensions, 
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already a gauge field configuration which behaves as 1/r in the vicinity of a 
point singularity will induce an action diverging as / dr r 3 (l/r 4 ); a divergence 
of the same degree is induced by a l/y/r line singularity or a ln(r) sheet 
singularity in four space-time dimensions. 



3.5.3 Alternative vortex gauges 

The remarks above illustrate that the continuum maximal center gauge in 



the form derived in section |37| in a sense constitutes excessive rigor. In a dy- 
namical calculation, one is forced to introduce a regularization of divergences 
in the ultraviolet; this regularization in particular also smears out the thin 
vortex configurations a(<9£), used in maximal center gauge fixing, into thick 
vortices with a thickness related to the ultraviolet cutoff. In view of this reg- 
ularization, dictated by the dynamics, it is not mandated to consider at an 
intermediate stage the behavior of gauge field configurations as the cutoff is 
removed. Instead, the cutoff must be kept finite throughout the calculation, 
and only the continuum limit of physical observables, as encoded in their 
renormalization group behavior, is relevant. 

This suggests an alternative formulation of the gauge fixing procedure in 
which one approximates a given gauge configuration as well as possible by 
thick vortices. Specifically, one can generalize the gauge condition ( j36|) and 
the remaining maximization problem (|88|) by allowing the approximating 
configurations a(<9£) to comprise thick vortices. Of course, such a general- 
ization clouds the original idea of the maximal center gauge, namely keeping 
the center part of the gauge freedom on the lattice unfixed; instead, it op- 
erates purely with continuum, i.e. coset configurations. Therefore, such 
gauges should be more appropriately called maximal vortex gauges rather 
than center gauges. 

In such vortex gauges, one has an additional freedom in the choice of the 
thickness of the approximating vortex configurations a(<9£). On the one 
hand, one may choose the thickness to be related to the ultraviolet cutoff 
with which one regulates the theory; for example, the thickness may simply 
be the lattice spacing. This form would be entirely equivalent to the maximal 
center gauge in the lattice theory, since, at the level of the regularized theory, 
there is no distinction between such a thick vortex and a truly thin one. On 
the other hand, one may choose the thickness to be a fixed physical quantity^ 

10 Some phenomenological implications of a physical vortex thickness were discussed in 

MM- 
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(which can be varied to optimally match the vortex profile preferred by the 
Yang-Mills dynamics). In this form, the thick vortex profile can also take 
over the role of the gauge fixing function g. For example, as evidenced in 
the examples given further above, the gauge fixing function g can act as a 
regulator for the divergences introduced into the gauge fixing functional by 
the singularities of the thin vortices. In a vortex gauge with thick vortices, 
this is already achieved by the vortex thickness. In many respects, the vortex 
profile in the vortex gauge has a similar effect as the gauge fixing function g 
in the maximal center gauge. 

As an aside, note that in a consistently regularized theory, one should also 
use regularized versions of the ideal vortex configurations A, cf. (p3|), which 
one simply obtains by using the ultraviolet cutoff of the theory to smooth 
out the 5-functions in ([33|), cf. eq. ( |133| ). 



Lastly, a certain aspect of continuum thin (or, equivalently, ideal) vortices 
has not yet been discussed in detail, namely the possibility that there may 
exist distinct thin vortex fluxes which contribute an identical phase to any 



Wilson loop they link. An example of this was given already in section ^72 
in the SU{2) thin vortex configuration, eq. (§), 

a„, = ±-T 3 , a r = (95) 
r 

both choices of sign lead to the same effect on a linked Wilson loop. Both 
choices correspond to a Cartan chromomagnetic flux described by E 3 = 2tt, 
cf. eq. (0), but they differ in the direction, or orientation, of the flux, which 
is inverted when the sign is reversed. 

On the other hand, on a Z(2) lattice, one cannot distinguish between these 
two possibilities, since one formulates the theory directly in terms of group 
elements. On a Z(2) lattice, there is only one type of center vortex flux, 
defined by plaquettes taking the value (—1). Information about the orienta- 
tion of the flux is lost during lattice center projection. For the purpose of 
evaluating Wilson loops, this is immaterial; however, there may be other ob- 
servables which are sensitive to the orientation of vortex flux. An important 
example of this, namely the Pontryagin index, will be discussed in detail in 
section |^. 

Thus, to obtain a comprehensive description of infrared phenomena in Yang- 
Mills theory, it is necessary to project onto a slightly more general class of 
infrared effective degrees of freedom than the essentially unoriented center 
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vortices defined by lattice center projection. Namely, vortex surfaces should 
additionally be associated with an orientation, and should in general consist 
of patches of differing orientation. As will be shown in section |5| in detail, the 
edges of these patches can be associated with Abelian magnetic monopole 
trajectories. The continuum center vortices introduced in the past sections do 
already include information about the orientation of magnetic flux. E.g., in 
the discussion after eq. (|76|) , when maximization over all hypersurfaces S is 
called for, different orientations of £ can be used. Likewise, in the discussion 
after eq. (0), different orientations of <9£, i.e. directions of flux in a, are 
distinguishable. However, vortex surface patches of alternating orientation, 
bounded by Abelian monopole currents, have hitherto not been allowed. On 
the contrary, the magnetic flux carried by thin vortex configurations a has 
been assumed to be continuous throughout. 

The purpose of the the present discussion is merely to round out the treat- 
ment of maximal center gauge fixing and center projection by pointing out 
how it can easily be generalized to yield vortex surfaces made up of ori- 
ented patches. More detailed properties of such surfaces, including how they 
determine the Pontryagin index, are discussed in section |5|. 

One way of introducing patches of different orientation on the vortex surfaces 
is to include information on Abelian monopole degrees of freedom into the 
gauge fixing and projection procedure. As mentioned above, the trajecto- 
ries of such monopoles on vortex surfaces define the edges of the oriented 
patches. Such a procedure has already been defined and implemented in 
the lattice formulation, namely the "Indirect Maximal Center Gauge" intro- 



duced in [10] and further investigated in [18|. This procedure is defined by 
initially transforming a given lattice Yang-Mills configuration to the maxi- 
mal Abelian gauge and performing Abelian projection, which allows to ex- 
tract Abelian monopole trajectories. In a subsequent step, one transforms 
the residual Abelian lattice configuration to the maximal center gauge and 
performs center projection, allowing the extraction of vortex surfaces. This 
two-step procedure, including two truncations of the theory by projection, is 
easily translated into the continuum formulation; the maximal Abelian gauge 
corresponds to the gauge fixing condition 

[dr + AM,Af>]=0 (96) 

where Ajf* denotes the color neutral (Cartan) components of the gauge field, 
and A^ denotes the color charged components. Fixing to this gauge induces 
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Abelian monopoles in the resulting field configurations. After projecting the 
gauge-fixed configurations onto their color neutral part (Abelian projection), 
the maximal center gauge fixing procedure, as described in the past sections, 
can be implemented analogously, and the associated vortices extracted. 

Note that the construction does not guarantee a priori that the monopoles 
are located on the vortex surfaces. However, monopoles have been found 
empirically to lie on the center vortices extracted via the indirect maximal 
center gauge in lattice experiments |Tj|. 

On the other hand, it is not strictly necessary to resort to such a two- 
step procedure, one of the steps being geared to extract Abelian monopoles, 
and the other to extract center vortices. The coset gauge transformations 
G in the maximal center gauge fixing functional fl82"D in particular include 
gauge transformations which generate Dirac strings^. Dirac strings, which 
describe open two-dimensional world-sheets in four-dimensional space-time, 
carry quantized magnetic flux such that they do not influence any Wilson 
loop, in accordance with the fact that they are unobservable pure gauge 
artefacts. E.g., in the case of SU(2) color, a Dirac string carries twice the 
magnetic flux of a thin vortex. Before discussing the interplay between vor- 
tices and Dirac strings, consider the following example for illustration. The 
SU (2) pure gauge configuration 

cos(£/2) sin(0/2)e-^ \ 
-sin(0/2)e^ cos(#/2) j K > 

where (r, 9, ip) denotes the usual three-dimensional spherical coordinates, and 
the configuration may be considered constant in the further space-time co- 
ordinate, has the explicit components 

) (98) 

nd irac i ^(Q/.J MO/2) COs(0/2) e -^ \ 

% - rsine 3111 ^/^ cos(0/2)e* -sin(0/2) J 1 ] 

11 Note that the description in terms of Dirac strings, used in the following for conve- 
nience, is equivalent to a description in terms of Abelian monopoles. The two objects are 
strictly coupled. Monopole trajectories define the edges of Dirac string world-sheets; Dirac 
string world-sheets span areas circumscribed by monopole trajectories, where the specific 
area associated with a monopole trajectory can be modified by gauge transformations. 



a dirac = GdQ ] ? q = 
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The full nonabelian field strength of this configuration can easily be verified 
to vanish p8 |, except on the negative z-axis0, i.e. for 9 = n. There, the 



Dirac string is located. For 9 — ► ir, a™ mc can be approximated by 



a dirac ^ 1 . = 2a thin vortex 

V r sin 9 v 

and the flux through an infinitesimal Wilson loop encircling the Dirac string 
therefore is 



2t 



$= / d<prsm9a d J rac = -4ttT 3 (101) 
Jo v 

leading, upon exponentiation, to the trivial value W — 1 for the Wilson loop. 
Furthermore, by considering the Abelian projection a dtrac of a dirac , 



—dirac 2 sin (0/2) —dirac —dirac r\ f 1 no\ 

a w — J-3 , 0-9 — a r — U 5 l iU/: J 

v rsm9 



one obtains in the associated Abelian magnetic field 

1 



B = d x a dtrac = -^e r T 3 (103) 



a Dirac magnetic monopole located at the boundary of the Dirac string. 
Note that in the full field strength, this field is completely canceled by the 



nonabelian part as mentioned above. Thus, magnetic monopoles which 
become visible in Abelian-projected configurations can be used as an alter- 
native way of detecting, or describing, the Dirac strings induced by coset 
gauge fixing. 

The above properties imply that a Dirac string world-sheet which coincides 
with a (oppositely oriented) patch of a thin 577(2) vortex world-sheet simply 
reverses the magnetic flux associated with that patch, leaving all Wilson loops 
invariant. Nevertheless, in this way, though representing a pure gauge on its 
own, a Dirac string may change what information contained in a configuration 
is kept during continuum center (or vortex) projection, i.e. it can be used to 
define a different truncation of the theory. The physical content of the full 
configurations of course is not influenced by the Dirac string. 

It should be noted that for higher SU (N) color groups, the superposition of 
a Dirac string on a vortex changes not only the sign (i.e. the orientation), 
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The action of this pure gauge configuration of course must vanish i n sp ite of the 



singularity on the negative z-axis, cf. the pertinent comments in Appendix A.l 
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but in general also the type of vortex flux, i.e. its direction in color space. 
E.g. for SU(3), one in general additionally exchanges the center flux labels 
(k = 1) <-> (k = 2) in eq. fl4Cf) . This happens precisely in such a way as to 
leave all Wilson loops invariant, as it should be. 

Thus, due to the possibility of vortex surfaces being partially covered by Dirac 
string world-sheets generated by coset transformations G in fl32|), vortex sur- 
faces made up of patches of different orientation are automatically generated 
during continuum center gauge fixing. Vortex projection must merely be gen- 
eralized to include this additional information about the orientation of vortex 
magnetic flux resulting from the Dirac world-sheets generated by coset gauge 
fixing. 

Finally, it should be noted that thick analogues of vortex surfaces consisting 



of patches of different orientation have been constructed in [2l^],||3~0|. Using 
such generalized objects as approximating configurations a in the gauge fixing 
condition fl82|) provides another, alternative, way of extracting information 
about vortex orientation from gauge field configurations. 



4 Sketch of an Effective Vortex Theory 

The center dominance observed for the string tension in lattice Monte Carlo 
calculations using the maximal center gauge |T?|,pT|,[|I0|,[|n],[]r2"| supports 
the notion that center vortices are the infrared degrees of freedom relevant for 
confinement in Yang-Mills theory. On a formal level, center or vortex dom- 
inance implies that the confinement properties can be adequately described 
by an effective vortex theory which results by appropriately integrating out 
all other Yang-Mills degrees of freedom. More specifically, in the abovemen- 
tioned lattice calculations, center dominance is obtained using the full lattice 
Yang-Mills dynamics, i.e. the field configurations are sampled according to 
the full Yang-Mills action. Center projection is only performed in the observ- 
able, i.e. the Wilson loop is evaluated using center-projected configurations 
after maximal center gauge fixing. Thus, in principle, one could evaluate the 
expectation value of the Wilson loop using a Z(N) effective action obtained 
after integrating out all coset SU(N)/Z(N) degrees of freedom. In lattice 
calculations, this separation of the functional integrations is impractical, and 
the Z(N) effective action is unknown. 

Adopting a continuum description does not greatly alleviate this problem, 
and therefore the following considerations concerning the effective theory of 
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vortices must remain mostly on a formal level. The intention is to sketch 
what types of effects may arise in the effective vortex theory particularly in 
the long-wavelength domain. 

To this end, it is necessary to implement a maximal center gauge in the 
Yang-Mills functional integral. This gauge fixing can be decomposed into two 
steps, as shown in section |3.4j , namely coset gauge fixing in an arbitrary fixed 
thin vortex configuration a and optimization of the thin vortex configuration. 
Coset gauge fixing can be achieved in the standard way via the Fadeev-Popov 
procedure, i.e. the gauge fixed partition function takes the form 

Z = J [DA] A F [A, a] exp (-S YM [A] - pF 2 [A, a]) (104) 

where F[A, a] denotes the gauge fixing condition ( jS6f) , Ap[A, a] is the stan- 
dard Fadeev-Popov determinant corresponding to this gauge fixing condition, 
and p is a gauge fixing parameter which should be sent to infinity in order 
to implement the coset gauge fixing condition ([SB]) exactly. Note that at this 
point a is an arbitrary, but fixed, external (thin vortex) field configuration. 
In order to optimize the choice of thin vortex field a, it is useful to further 
introduce into the functional integral the following representation of unity^ 

1 = A[A] J[DE\exp(-S G p[A,a(dE)]) (105) 
S GF = -q J d D x g (tr(A - a) 2 ) (106) 

cf. eq. flgE|). Here, the integration runs over all (D — l)-dimensional hyper- 
surfaces £ in D space-time dimensions, and the specific realization a(<9£) of 
a thin vortex on <9£ can e.g. be given by (|6~2|) . Furthermore, the gauge fixing 
parameter q should also be sent to infinity in order for e~ s ° F to be peaked 
at the optimal thin vortex configuration. 

Note that ( |105| ), in order to be properly defined, at the very least requires 
some further gauge fixing due to the fact that the integrand depends only on 
the boundaries <9£ of the hypersurfaces S being integrated over. Note that 
it is by no means clear that the integral over hypersurfaces £ can simply be 
replaced by an integral over closed surfaces S; there may exist closed surfaces 
S which cannot be represented as S = <9£ in terms of some hypersurface S. 

13 Note that one could also forego separating coset gauge fixing and vortex optimization, 
and instead use the full continuum maximal center gauge fixing functional (^2|) to introduce 
a representation of unity into the Yang-Mills functional integral analogously to ( |105[ ) . 
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A simple example is one of the coordinate planes, e.g. the 1-2-plane, in 
a four- dimensional space-time hypercube endowed with periodic boundary 
conditions, i.e. a torus. Such a plane constitutes a closed surface due to the 
periodic boundary conditions, but cannot be represented as the boundary 
of any three-volume S. In general, the number of topologically inequivalent 
closed two-dimensional surfaces which are not boundaries of three- volumes is 
given by the dimension of the 2nd homology group ^(M) of the space-time 
manifold M under consideration; this corresponds to the 2nd Betti number 
b 2 = dimH 2 (M). On the four-dimensional torus T 4 , one has &2(^ 4 ) — 6, 
corresponding precisely to the six independent planes which are closed by 
the periodic boundary conditions. 

Formally, one can e.g. modify ( |105| ) by writing 

1 = A[A] J [DS] J [DE]5[E - E min [S}} exp (S GF [A, a{S)\) (107) 

where S m j n [S'] denotes the hypersurface of minimal volume which has S as 

its boundary, <9£ min = S. If no hypersurface E exists such that <9£ = S, then 
that S does not contribute in (|107|) . Note that this modification, leaving from 
the original integration over hypersurfaces £ only an integration over closed 
(D — 2)-dimensional surfaces S in D space-time dimensions, corresponds 
precisely to a fixing of the center part of the gauge freedom on the lattice. 
Note also that this further gauge fixing in general leads to further measure 
factors, cf. the modified measure A[A] in (|107|) as opposed to A [A] in ( |105| ). 

Inserting (|107Q into ( |104j) , and assuming functional integrations can be in- 
terchanged, one arrives at 



Z = f[DS}e- s ^f [s] (108) 



J[DZ]5[X-Z min [S]] (109) 

'[DA] A[A}A F [A, a } e SYM[A]-pF2[A,a]-Se F {A,a] 



where a = a(S), as specified e.g. by (|62|). The effective action given by 
( |109| ) describes a theory of (D — 2)-dimensional closed vortex surfaces 5* in 
D space-time dimensions. In the most general case, e.g. for the purpose of 
evaluating the vortex-projected Pontryagin index, cf. next section, one will 
supplement the surfaces S with information about Dirac strings contained 
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in the coset gauge-fixed configurations A; this leads to patches of different 
orientation on S as indicated in section |3.5.3| . Formally, the surfaces S being 



integrated over in (|108|) will consist of oriented patches, and the integration 
over A in (|109| ) for given S will in turn be restricted to all A containing the 
corresponding Dirac string configuration. 

Of course, the previous manipulations are at this stage purely formal. Even 
to obtain the measure ALA] in ( |109| ), one has to solve a string theory in an 
arbitrary gauge background A in four space-time dimensions, namely one 



has to carry out the integral over S in (|107|) . However, in lower dimensions, 
an evaluation of ( [LU7| ) may be feasible. E.g., in three space-time dimensions, 
the S denote closed lines, and the integration over S may be converted to a 
field theory by interpreting these lines as world-lines of particles. 

In order to nevertheless gain some insight into the structure of the vortex 
effective action at least in the infrared limit, assume the validity of the state- 
ment that the center projection vortices in Yang-Mills theory give a rough 
localization of thick vortices present in the full Yang-Mills configurations. 
As discussed in section [2.2|, this statement is corroborated by empirical find- 



ings in lattice experiments |T0[ and can be understood on the basis of the 
continuum formulation of the maximal center gauge developed in section [3[ 
However, it is a statement which can at most be considered valid on length 
scales coarser than the typical thickness R of the aforementioned vortices in 
the full configurations preferred by the Yang-Mills dynamics. At this point, 
R is an undetermined scale; it should ultimately be extracted by minimizing 
the effective thick vortex action, as indicated further below. 

The above statement implies that the distribution of thin projection vortex 
configurations in the infrared regime is simply determined by the distribu- 
tion of thick vortices in Yang-Mills theory. In other words, the infrared 
effective action of projection vortices can be obtained by calculating the 
standard Yang-Mills effective action of thick vortex configurations, using a 



thick vortex background gauge as already briefly discussed in section [3.5.3 
The ultraviolet limit of validity of this effective action is given by the cutoff 
A = l/R. 

Furthermore, in such a long-wavelength approximation scheme, it is consis- 
tent to evaluate the effective action in a gradient expansion, i.e. in powers 
of Rdi, where di = djdoi is the derivative with respect to the parameters of 
the vortex surface. Below, this effective action will be given in the classical 
limit, i.e. to zeroth order in a loop expansion. This will already exhibit the 
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relevant terms, with higher orders in the loop expansion only renormalizing 
the coefficients. The calculation is carried out explicitly in Appendix |A.1| us- 



ing, for convenience, ultraviolet-regulated ideal vortex configurations. More 
precisely, the 5-function in the defining expression (^) for the ideal vortex 
A(k, S, x) is regularized by means of the cutoff A = 1/R, cf. eq. (|133|) . Note 



that such an object represents a bona fide thick vortex; the special character 
of ideal vortices as singular continuum representations of excitations on a 
Z(N) lattice is blurred when the ideal vortex singularities are regularized. 
After regularization, there is no qualitative distinction between an initially 
thin or an initially ideal vortex; these can all be viewed as continuum coset 
configurations, namely, thick vortices. 

Inserting such a thickened ideal vortex A into the classical Yang-Mills action, 

S[A\ = ~ ftr(F[A)F[A)) (110) 
one obtains the gradient expansion (see Appendix |A.1| ) 



Sym[A] =a [ d 2 a^ + (3 f d 2 a^K^g ab g^ 

JdT, JdT, 

+ 7 1 d 2 a^gg ab d a d b ln^g . (Ill) 

Here, 

g = detg ab , g ah = daX^d^ (112) 

denotes the metric on the vortex sheet and K^ b are the extrinsic curvature 
coefficients which are defined in Appendix [A.l| , eq. (|138|) . Furthermore, 



the coefficients in eq. ( |111 ) are given for SU(2) vortices by (using trE" 2 
4(27r) 2 , cf. eq. (H) ) 



a =Jp- "=-^- ^ii- (113) 



One thus obtains in leading and next-to-leading order the Nambu-Goto action 
and curvature terms for the vortex surfaces^. The coefficients ( |113| ) will be 
renormalized through quantum effects embodied by higher orders in the loop 
expansion. In particular, renormalization will generate a physical scale in the 
coefficients in the standard fashion, which will allow to relate the thickness 
R to physically measurable quantities. This will make it possible to self- 
consistently determine the preferred physical vortex thickness R through a 



J Thc second term on the r.h.s. of eq. (|1 1 1|) has been discussed in ]31|,|3 



43 



minimization of the effective action for given vortex magnetic flux; similar 
considerations were already carried out in the framework of the Copenhagen 
vacuum 0. Such more detailed approximations, e.g. an evaluation of the 
one-loop corrections, are deferred to future work. 

It should however be emphasized that, though the explicit form of the coeffi- 
cients of the various terms in the effective vortex action certainly depends on 
the approximation involved, the general structure of the gradient expansion 
of the vortex effective action is uniquely determined by Poincare, gauge, and 
reparametrization invariance. In a forthcoming paper, the authors will use 
an action reminiscent of (|1 1 1|) with adjustable coefficients as a starting point 



for a phenomenological vortex model of infrared Yang-Mills physics. 



5 Topology of Center Vortices 

In the previous section, a rough sketch of a reduction of Yang-Mills theory to 
an effective vortex theory was given, by implementing a continuum version 
of the maximal center gauge, and by subsequently integrating out the Yang- 
Mills fields. Center-projected (or, more appropriately, vortex-projected) ob- 
servables can in principle now be evaluated directly in terms of the vortex 
theory, which contains exclusively (D — 2)-dimensional vortex surface degrees 
of freedom in .D- dimensional space-time. 

There is an essential difference between the center vortices on the lattice and 
in the continuum theory. After center projection on the lattice, the direction 
of the magnetic flux of the vortex is lost, while in the continuum theory, the 
center vortices are given by oriented (patches of) surfaces, where the flux 
direction is defined by the orientation. This was already discussed to some 
extent in section |3.5.3| . The orientation of the vortices, however, is crucial 
for their topological properties, as will be shown below by evaluating the 
Pontryagin index 

M\ = 3^2 / d ^K F U = / d 4 xtvF^F, u (114) 

using the ideal vortex representation (P5|). 

It should be emphasized that the treatment in this section, the goal of which 
is to define the vortex-projected Pontryagin index, in some ways reverses the 
logic of the past sections. Whereas the continuum maximal center gauge and 
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also the vortex action were defined via the corresponding lattice expressions 
in the limit of vanishing lattice spacing, the Pontryagin index calls for a 
different treatment. It is initially defined for smooth continuum gauge field 
configurations, and useful lattice discretizations of the Pontryagin index are 
notoriously cumbersome. In fact, the maximal center gauge exacerbates the 
problem, since all lattice definitions at some point must assume the links of 
a lattice configuration to be reasonably close to unity. On the other hand, 
the maximal center gauge precisely transforms configurations into a form 
where their links strongly deviate from unity; namely, it introduces ideal 
vortices even into initially smooth configurations. Therefore, to properly 
define the vortex-projected Pontryagin index, it is necessary to start from 
the continuum expression, understanding the Pontryagin index for thin (or, 
gauge-equivalently, ideal) vortices to correspond to the thin limit of the Pon- 
tryagin index of thickened vortices, for which it is well-defined. This limit is 
smooth, as is demonstrated in an example in Appendix [C|; namely, the thick 
vortex profile function manifestly cancels. Note that this concept of vortex 
projection is not directly connected with lattice center projection anymore, 
since it is based on the thin limit of thickened (i.e. coset) configurations. In 
the following, the Pontryagin index will be evaluated for (thickened) ideal 
vortex configurations. The aforementioned example in Appendix |C| on the 
other hand treats a particular configuration of explicitly thick vortices, con- 
taining also the limiting case of thin vortices. 

The field strength of an ideal center vortex A is given by 

F^(k,x,S = dE)=F ia ,[A] = d li A v (k,E,x)-d u Ar(k,E,x) (115) 

Straightforward evaluation (see Appendix |A.1| ) yields for the field strength^ 
on a vortex surface S 

F^k, x, S) = E{k) f d 2 a^ v S 4 (x - x(a)) . (116) 

One can decompose a general vortex configuration into components of differ- 
ent center flux labeled by the integer k, where k = 1, . . . , iV — 1 for SU(N) 
color, 

S = \JS k , F/wfaS) = '£j 7 flu/ (x,Sk) (H7) 
k k 

Inserting this, and eq. (|116p , into (|114 ), one finds for the Pontryagin index 



v[A] = --^ J (Sk, S k ,)tr(E(k)E(k')) (118) 

° 7F k,k' 
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In the context of superfluid Helium, this quantity is referred to as the vorticity tensor. 
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where 

I{Si, S 2 ) = \ [ d 2 a,„ [ d 2 a'^S 4 (x(a) - x(a')) (119) 

is the (oriented) intersection number between two two-dimensional surfaces 
Si, S 2 in D — 4 space-time dimensions. Note that the intersection number 
I (Si, S 2 ) is also defined for open surfaces Si,S 2 , and that, in D = 4 dimen- 
sions, /(Si, 62) = I(S2,Si). Also, if Si and S 2 intersect in a single point, 
then I (Si, S 2 ) = ±1, where the sign depends on the relative orientation of Si 
and S 2 . Note furthermore that, in the case k = k' in ( |118| ), I(Sk, Sk) counts 
each isolated self-intersection point of Sk twice. This is easily seen by further 
decomposing S fe into two components Si and S| such that neither of the two 
components self-intersect at the point in question. Then, 

I(S k , S k ) = I(Sl Si) + I(S 2 k , Si) + 2I(Sl, S 2 k ) (120) 

where the final term represents the contribution from the self-intersection 
point of Sk under consideration. 

It should be remarked that this discussion of the self-intersection number 
I(S, S) includes only the contributions from transversal intersection points, 
i.e. points where x^(a) = x^(a') with a 7^ a', but leaves out the contribution 
from the coincidence points x^(a) = x^(a') with a = a'. A more thorough 
discussion of the self-intersection number I(S, S), which is given in Appendix 
[B], does, however, not alter the conclusions on the Pontryagin index reached 
here. 

Thus, the Pontryagin index of a vortex surface configuration is built up from 
its self-intersection points, each of which enters with a contribution ±2 • 
tr(E(k)E(k')) /8tt 2 into the Pontryagin index according to ( |118| ). The integers 
k, k! label the center flux carried by the two patches of the vortex surface 
configuration intersecting at the point in question; the sign depends on the 
relative orientation of the patches. In view of eq. (|38|), this implies specifically 
for SU(2) color a contribution ±1/2 from each self-intersection point; on the 
other hand, for SU(3) color, each self-intersection point contributes ±1/3 or 
±2/3 to the Pontryagin index, cf. eqs. ( |41~|) and fl42|) . This is consistent with 



the dependence on the number of colors noted by J. M. Cornwall in p9 



In view of this result, the reader might worry at this point about the integer- 
valuedness of the Pontryagin index on R A (silently compactified to S A ). This 
property cannot go lost during vortex projection, for the very simple reason 
that vortices can be explicitly represented as gauge field configurations and 
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therefore inherit all properties known to be true for generic gauge fields. 
Nevertheless, it is instructive to remark how the integer-valuedness comes 
about in more detail in a simple case, namely SU(2) vortices on i? 4 . It is due 
to the fact well-known in topology^ that the number of intersection points of 
two closed two-dimensional surfaces in R A is even. This also implies that the 
number of self-intersection points of a closed surface is even, because the self- 
intersection number is defined by simply intersecting the surface with another 
surface infinitesimally displaced from it (i.e. one considers a framing of the 
surface, cf. Appendix [B]). Thus, while each self-intersection point of an SU(2) 
vortex surface configuration gives a contribution ±1/2 to the Pontryagin 
index, the number of such contributions is even, and the Pontryagin index 
integer-valued. 

In the case of higher SU(N) gauge groups, the argument becomes more 
complicated, since the surfaces may branch and also due to the fact that 
superimposed Dirac strings in general also modify the type of vortex flux, 
i.e. its direction in color space. Further complications arise on different 
space-time manifolds. For instance, on a four-cube with periodic boundary 
conditions, i.e. a torus, one might think of generating a single intersection 
point between two surfaces by choosing e.g. the 1-2 plane and the 3-4 plane, 
which are closed due to the periodic boundary conditions. However, this 
example is invalid, since vortex surfaces must be representable as surfaces of 
three-volumes, which is not the case for the aforementioned planes. 

The final point to be discussed is the fact, also well-known in topology, that 
the self-intersection number of closed, globally oriented two-dimensional sur- 
faces in R 4 not only is a multiple of four (remember that each intersection 
point is counted doubly, and there are an even number of intersection points 
for closed surfaces), but actually vanishes. This implies that the Pontryagin 
index vanishes for globally oriented vortex surfaces; conversely, therefore, 
non-orientedness of the surfaces is crucial for generating a non-vanishing 
topological winding number. To illustrate this, consider the following three- 
dimensional slice of two intersecting SU{2) vortex surfaces in D = 4 dimen- 
sions, cf. Fig. H Let the one vortex be located entirely within the three- 
dimensional slice of space-time under consideration; it is therefore visible as 
a closed surface, namely the sphere S in Fig. |3|. On the other hand, let the 
other vortex extend into the space-time dimension not displayed in Fig. ^ 
it is then visible as a closed loop C after slicing. Let C intersect S at two 

16 The authors thank B. Leeb for a discussion on this point. 
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points, chosen in Fig. [5] as the poles of S. 



c 



c 




s 




(a) 



(b) 



Figure 3: Three-dimensional slice of two intersecting SU(2) vortex surfaces 
in four space-time dimensions (see text). The thin line around the equator 
of the sphere S in (a) is merely to guide the eye in identifying S as a sphere. 
At the intersection points of S and C, contributions to the Pontryagin index 
of modulus 1/2 arise. However, in the case of globally oriented surfaces (a), 
these contributions cancel. On the other hand, in the case that vortex sur- 
faces are not globally oriented, but consist of patches of different orientation 
(b), the contributions may add, giving rise to a nonvanishing Pontryagin in- 
dex. Boundaries of such patches are tagged by Abelian monopoles (thick line 
around the equator in (b)). 

If S and C are both globally oriented, as in Fig. [3] (a), their intersection 
number (in D = 3) 



vanishes, since the two intersection points occur with opposite relative ori- 
entation between the intersecting manifolds. Thus, a non-zero intersection 
number requires at least one of the two intersecting closed manifolds to be 
not globally oriented. Globally non-oriented hypersurfaces consist neverthe- 
less of oriented patches. Assume for example the sphere S in Fig. |3| (b) to 




(121) 
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consist of two hemispheres of opposite orientation. E.g., the southern hemi- 
sphere may be covered by a Dirac string world-sheet, simply reversing (in 
the case of 577(2) color) the orientation associated with that hemisphere. 
The boundary of a Dirac string world-sheet is complemented by a monopole 
in the Abelian-projected configuration, cf. eq. (|103|) . Thus, in Fig. § (b), 
one can associate an Abelian monopole current with the equator (thick line). 
Further below, it will be shown more generally how boundaries of oriented 
patches on vortex surfaces imply monopole currents. 

Now, in the presence of the monopole current at the equator in Fig. [| (b), 
implying a reversal in orientation of, say, the southern hemisphere, the con- 
tributions to the Pontryagin index from the two intersection points no longer 
cancel, but add coherently, giving a Pontryagin index of modulus \v\ = 1. 
Thus, the depicted configuration represents an analogue of an instanton in the 
sense that it will generate a zero mode in the Dirac operator according to the 
index theorem. The authors speculate that standard smoothing procedures 
used on the lattice to detect instantons would transmute this configuration 
into an instanton, since smoothing by construction alters configurations such 
as to lower the Yang-Mills action. After smoothing, it is impossible to dis- 
cern what type of excitation precisely carried Pontryagin index in the original 
configuration. Note that this picture of the generation of topological winding 
number also correlates well with the observation |33|],[]34j that instantons are 
encircled by monopole loops in the maximal Abelian gauge. It is precisely 
the monopole loop which tags the change in orientation in the vortex surface 
(the sphere 5 in Fig. ^ (b)) between the two intersection points carrying the 
topological density, causing their contributions to the Pontryagin index to 
add up and generate the "pre-instanton" . 

To see the emergence of the magnetic monopole loops at the boundaries of 
vortex surface patches in more detail, consider an isolated patch; this now 
represents an open two-dimensional surface sheet 5 with boundary dS = C. 
To this surface 5, one can relate a magnetic flux sheet by associating with it 
the field strength (see above) 

F^x, S)=E f d 2 a^5\x - x{a)) . (122) 
J s 

For a closed surface S = <9£, this expression agrees with the one for the 
ideal vortex „4(E,x) given in (|116|). By Stokes' theorem, one finds that a 



conserved (Abelian) magnetic current 



ff(x, C) = d v T, v = El _ dy,5\x - y) , d^(x, C) = (123) 

•J OS — c 
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must flow at the boundary dS — C of the magnetic sheet 5*. Note that the 
direction of the magnetic current on dS = C defines an orientation of the 
sheet S and vice versa. This monopole current is generated by a magnetic 
point charge q m 

Eq m = f da^(x, C) = EI(V, C) (124) 

with V denoting the three-volume dual to the loop C. In view of exp(— E) = 
Z, eq. ( |124| ) yields for q m only a fraction of the elementary magnetic charge 
dictated by the Dirac quantization condition. This is as it should be, since the 
vortex sheet S only carries a fraction of the flux which a Dirac string carries. 
Only when glueing the open vortex patches with their boundary monopole 
loops back together, should one take care to obtain current loops of proper 
magnetic point charges satisfying the Dirac quantization condition. In this 
way, non-oriented closed magnetic vortex sheets are constructed, consisting 
of oriented surface pieces joined by magnetic monopole current loops^]. Such 
a non-oriented closed magnetic sheet defines a center vortex which still gives 
the same contribution to a Wilson loop as the corresponding oriented vortex 



(in the absence of a monopole loop), cf. eq. (43). Thus, for the confinement 
properties, measured by the Wilson loop, the orientation of a vortex sheet, 
and hence the magnetic monopole currents, are irrelevant. The Abelian 
magnetic monopole currents are however necessary in order to generate a 
non-vanishing Pontryagin index for vortex configurations. It is in fact easy 
to see that the Pontryagin index of the vortices can be determined from the 
Abelian magnetic monopoles alone. Indeed, inserting ( |115| ) into ( |114| ), after 



a partial integration the Pontryagin index of an ideal vortex sheet becomes 
v = / d fl tr(A l/ T fll/ ) + ^ [ trA<V^ ( 125 ) 

07T Z JM o7T z JM 

= ~ I ^r(XF ) + _L/ tr A„j?{x,C) (126) 

o7T JdM on 2 JM 

where Gaufi' law was used in the first term and the definition of the mag- 
netic current in the second term. The second term obviously vanishes in 
the absence of magnetic monopoles. Furthermore, the surface term van- 
ishes unless the integrand contains singularities which give rise to internal 
surfaces wrapping the singularities. Such internal surfaces precisely arise in 



the presence of magnetic charges [14] which may be magnetic monopoles or 



17 If one glues together two patches such that the surface orientation does not change 
across the boundary, the monopole currents at the boundary precisely cancel. 
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more extended magnetic charge distributions such as line or surface charges. 
Thus, a non-zero Pontryagin index requires the existence of magnetic charges 
in Abelian-projected configurations. This is consistent with the findings in 

0- 



6 Outlook 

In this work, the continuum analogues of the maximal center gauge and 
center projection in lattice Yang-Mills theory were constructed. This shed 
new light on the meaning of the procedure on the lattice and led to a sketch of 
an effective vortex theory in the continuum. Also the manner in which center 
vortex configurations generate the Pontryagin index was clarified; the latter 
is built up from self-intersections of the vortex network, where it is crucial 
that the vortex surfaces be globally non-oriented. Since these developments 
have already been discussed at length in the corresponding sections above, it 
suffices here to mention some open issues which require further investigation. 



With regard to the maximal center gauge, it is necessary to conduct lattice 
experiments using alternative gauge fixing functions g, cf. eq. @, and to 
test whether the physical conclusions reached using the direct maximal center 
gauge remain unaffected. If it is true that maximal center gauges provide 
a rough localization of thick vortices present in full lattice configurations, 
then long-distance properties should be robust with respect to a variation 
of g; fluctuations of the center projection vortices within the thick vortex 
profile are of course gauge- dependent, as evidenced in section |272|. An issue 
which is of particular interest in this context is the occurrence of Gribov 
copies which do not correspond to the absolute maximum of the center gauge 
fixing functional. In fl9|| , the possibility is raised that the true gauge fixing 
image selected in particular by the direct maximal center gauge contains 
less vortices than found in previous studies, and no vortex confinement. As 
discussed in sections |3.5.1| and |3.5.2| , the direct maximal center gauge in the 
continuum limit in fact always selects a trivial gauge fixing image devoid of 
center vortices. The observations of [T9j may be a signature of the specific 
problems of this particular gauge. On the other hand, it is straightforward to 
formulate alternative gauge fixing functions g which prevent these problems. 
It would be interesting to repeat the considerations of |T]1 using such a gauge. 
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Concerning the issue of an effective vortex theory, it was already mentioned 
in section [| that the authors plan to report on their numerical study of a 
phenomenological random surface model for vortices in a forthcoming paper. 
On the other hand, from an analytical point of view, it would presumably be 
helpful to start with lower-dimensional models. E.g., it would be worthwhile 
to study whether the maximal center gauge allows a solution of Yang-Mills 
theory in 1 + 1 space-time dimensions; there, the effective vortex theory takes 
the form of a classical statistical mechanics of vortex points. 

Finally, it is necessary to investigate whether center vortices do in fact fully 
capture the physics of the Pontryagin index; i.e. one should test whether 
the vortex-projected Pontryagin index displays vortex dominance in analogy 
to the center dominance observed for the string tension. While it has been 
shown in this work how vortices generate a nontrivial topological winding 



number, and the converse experiment in [16[] has shown that an ensemble 
devoid of center vortices is concentrated in the trivial topological sector, the 
question of vortex dominance is more stringent. In this context, it is worth 
mentioning the analogous result from Abelian-projected theories: In the max- 
imal Abelian gauge, one finds monopole dominance in the string tension M, 
but not in the Pontryagin index; on the other hand, in the Polyakov gauge, 
the monopole dominance in the string tension is much less pronounced, but 
the Pontryagin index can be reconstructed exactly from the monopole con- 
tent of the gauge-fixed configurations alone |L4] . 



In connection with the issue of the Pontryagin index, of course also a more 
detailed study of chiral symmetry breaking by a vortex background is called 
for. Vortex configurations which give rise to a non-vanishing Pontryagin 
index must, by the index theorem, give rise to zero modes of the Dirac 
operator. It would e.g. be interesting to explicitly construct such zero modes 
for some relevant simple vortex configurations. In the context of instanton 
models, the zero modes induced by the instantons give in fact the dominant 
contribution to the chiral condensate. The question arises whether the chiral 



condensate displays vortex dominance (again, the converse experiment |L6 
yields zero condensate in the absence of vortices), and, if so, how it behaves 
at finite temperatures, in particular, at the deconfinement transition. 
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A Properties of center vortices 
A.l The ideal center vortex 

In the following, the field strength of the ideal vortex field A(k, S, x) 

is calculated. For this purpose, it is convenient to evaluate first its dual, 

fnv = 2 e »™/3Fai3[A) = e^uapdaApih, E, x) (127) 

Inserting here the explicit representation for Ap{k, S, x), cf. eq. (|33"D , one 
obtains 

^{k,dT.,x) = E{k)e^ VKX Jj D - l d x dl5 D {x-x{a)) (128) 
= E(k)^ UKX e Xafry J d^a^d^ix-xia)) (129) 
= E(k) Jj D - l o K , u dl5 D (x-x(a)) (130) 

Application of Stokes' theorem yields 

^(k, <9£, x) = E(k) f d D ~ 2 o^5 D {x - x(a)) . (131) 

Taking the dual of this equation yields the desired representation (|116|) . 

With the vortex field strength eq. ( |131| ) the Yang-Mills action of an ideal vor- 
tex in D = 4 space-time dimensions is given (up to an unimportant numerical 
factor) by 

S=f d 2 a, u f dV/(x( ff )-xK)). (132) 



.as ' Jas ^ 

Obviously this action is divergent and calls for regularization. Regularization 
will be achieved by chopping off the high frequencies using an ultraviolet 
cutoff A; this amounts to replacing the 5-function in equation ( |132|) by the 
regularized version 

5(x) = -^=e~^ 2 . (133) 

V Z7T 



Before proceeding, a comment is in order concerning the thin limit A — > oo. 
As long as A is arbitrarily large, but fixed with respect to the ultraviolet 
regulator of the theory, say the inverse lattice spacing, it is legitimate to 
evaluate the action using the continuum Yang-Mills expression (|132|) . How- 
ever, if the vortex is meant to be truly thin, i.e. if its entire magnetic flux is 
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concentrated on one lattice plaquette, then one must use the lattice plaquette 
expression for the action; otherwise, one loses track of the compact character 
of the gauge group. E.g., if one exactly superimposes two truly thin SU(2) 
ideal vortices on top of each other (this corresponds to a closed Dirac string), 
then naive application of the continuum Yang-Mills expression for the action 
would lead to a value diverging like four times the action of a single vortex. 
On the other hand, all plaquettes encircling this Dirac string take the value 
(+1), and its action therefore vanishes, as it shoulcQ. 

Keeping this caveat in mind, the regularized action can be evaluated in a 
gradient expansion. Note that similar calculations have been performed in 
[j35|, [|36| for a massive propagator. Since the action receives contributions 



only from x(a) ~ x(a'), the integrand ( |132| ) can be expanded in powers of 
z — a — a'. Defining s = and expanding (where d a = d/da a ), 

z a 

x(a) = x(s) + —d a x(s) H 

z a 

x(a') = x(s) d a x(s) 



2 

z a 1 z a z^ 

E jUi ,(cr) = E^(s) + — dalles) + -— —d a d b Y,^(s) + 



= - y + ~y3 ft ^E p (s) + • • ■ , (134) 



one obtains 



An / JdT, 



z a z b „ „ , . „ „ . . Z" z 



——d a Y, IJ , I/ ((j)d b Y, IJlu (a) + —— E IJ , u (a)d a d b E IJlu (a) 



(135) 



where only terms up to order z 2 in both the exponent and preexponent have 
been kept. Furthermore, the range of integration over z has been extended 
to ±oo, since the main contribution to the integral comes from the region 
zwO. Carrying out the Gaussian integration over z, one obtains 



dg a b y/g. 



(136) 



18 In the case of the Dirac string, these subtleties can be circumvented by dividing the 
space-time manifold into coordinate patches as in the Wu-Yang construction. 
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Here, the metric on the vortex sheet has been introduced, 
9db(p) = daX^dbX^a) , g = det(g ab ) . 



(137) 



Furthermore, a partial integration has been performed in the last term of 
equation ( |135| ). The expression obtained above can be cast into a more 
standard form by using the decomposition 

d a d b x,(a) = Y c ab d c x, + K&rf , (138) 



Furthermore, riA with A = 1,2, are unit vectors orthogonal to the tangent 



where T^ b are the affine connections and K^ b denotes the extrinsic curvature 

•more, raA with A 
vector of the vortex sheet 

n i d cX^{a) = , n*n 

Using 



\AB 



A = 1,2 . 



,-O-b cd 



9 ac -9 



(139) 
(140) 



one eventually obtains after straightforward evaluation 



S 



A 



2tt 



/ d 2 cjyfg 

JdT, 



1 

A2 



g ab g lk K*K£ k + ^g ah d a d b In ^ 



■ (141) 



The first term, which diverges for A — * oo, is the familiar Nambu-Goto term. 
The remaining two terms are finite for A — >• oo. The second term, containing 
the extrinsic curvature coefficients K^, has been discussed in fl31|,|j3^1. The 
higher-order terms in the gradient expansion all vanish for A — >• oo. 



A. 2 The thin center vortex 

In the following, it is shown that the thin center vortex as defined by eq. (|62]) 
for arbitrary vortex shapes 5* = <9£ indeed yields a center element for the 
Wilson loop when non-trivially linked to it. For this purpose, it is convenient 
to write the thin vortex as 

afl (k,S,x) = E(k)a^S,x) (142) 

a^S,x) = [ d D - 2 a^ K d x K D(x - x(a)) (143) 
J s 
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Using Stokes' theorem, one obtains 

f dxu,aJS,x) = / d 2 axu,dZaJS,x) (144) 

JC JM(C) 

= [ d 2 a x J d D -% K dld x K D(x - x(a)) (145) 

where C = dM. By using the relation between da^ and its dual da^, the 
last expression can be rewritten as 

dx^a^S, x) = -- ^ d 2 cr a/l3/ d D ' 2 a^e^ val pie^ Ka pd\d x K D{x - x(a)) 

(146) 

Using the properties of e a ^s, this expression can be converted to 

I dx^a^(S,x) = -- / d D ~ 2 a a p(x) [ d 2 a af3 {x)d 2 x D(x - x) (147) 
Jc 2 Jm(C) Js 

+ [ d D ~ 2 a K(3 { d 2 a aP d a d K D{x - x)(148) 

JA/(C) 

The last term is seen to vanish by Stokes' theorem, since dS = <9(<9E) = 0. 
The first term yields, upon using the definition of the Green's function, the 
intersection number between M(C) and S, 

I dxaaJS,x) = - [ d D ~ 2 a a p(x) [ d 2 a a p(x)6 D (x - x) (149) 

JC l JM(C) JS 

= I(M(C),S) (150) 
= L(C,S) (151) 

which equals the linking number between C and S, i.e. L(C,S). With the 
last relation, one indeed obtains a center element for the Wilson loop from 
the thin vortex (1621) , 



exp (- £ dx^k, S, x)) = exp(-E(k)L(C, S)) = Z(k) L{c ' s) (152) 

B Intersection number of two-dimensional 

sheets in D = 4 

This appendix summarizes essential properties of the self-intersection number 
of two-dimensional surfaces in D = 4 space-time dimensions; in particular, 
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the relation of the self-intersection number to topological invariants of the 
three-dimensional spaces located at the initial and final times is presented. 
Assume a certain parametrization of the vortex sheet = x M (o"i,cr 2 ) 

in D = 4. The self-intersection number I(S,S), cf. eq. ( |119| ), i.e., the 
intersection number of 5* with itself, receives contributions from "coincidence 
points" x m (<j) = x^a') with a = a' and from (transversal) "intersection 
points" x^{a) = x /t (cr') with o ^ a 1 , which will be denoted by Ii(S,S) and 
I 2 (S,S) respectively, 

I{S,S) = h{S,S) +I 2 (S,S) (153) 



Regularizing the delta function in eq. ( |119| ), cf. Appendix [A.l| , eq. 
introduces an ultraviolet cutoff A, which then allows I\(S, S) to be straight- 
forwardly evaluated; in the limit A — > oo, one finds Polyakov's intersection 
number^, 

h(S, S) = / f^g^d^d^ (154) 

167T Js y/g 

where 

9ij = diX ^((7)83X^(0) (155) 
denotes the metric on the surface S, g = det g^ and 

= tijdiX^^djXu^a) (156) 

is the tensor of the surface element on the vortex sheet. 

Let x M (cx(fc)) = x^(a'(k)), where a{k) 7^ cr'(fc), determine the coordinates of 
isolated (self-)intersection points P{k) of S. Expanding x M (<r), x^a') near 
a = cr(k), a' = <y'(k) in powers of a — cr(k), a' — cr'(k), the integrals in ( |119| ) 
over d 2 a, cPa' can be straightforwardly performed, yielding 

I 2 (S,S) = 2 ^ sign fc^(x(a(k)))£^(x(a(k)))) (157) 

k 

One observes that each intersection point x^aik)) = x^(a'(k)), with cr(fc) 7^ 
cr'(k), contributes doubly to the self-intersection number I 2 (S,S) in agree- 
ment with the findings in section |5|. 



19 



Ii(S, S) differs from Polyakov's original definition |u| by a factor of (—j 
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To illustrate the geometric meaning of the different contributions to the in- 
tersection number, it is useful to consider a fixed time slice of the four- 
dimensional universe. In this three-dimensional space, the vortex S repre- 
sents a closed loop C(t = fixed). For simplicity, choose t = a 2 such that this 
loop is parametrized by a±, which takes values in [0, 2ir[. 

Consider the self-linking number SL(C, h) of such a closed loop C(t) in D = 
3, which is defined by the linking number L(C, C) between the original loop 
C(t) and the fictitious loop C'(t) resulting from C(t) by an infinitesimal shift 

M°"i> *) + e ™(°"i> *) • ( 158 ) 

Here, 71(01, t) is a unit vector perpendicular to the tangent vector x (cri,t) = 
-^x(ai,t) of the loop C(i), 

h(aut)-x(ai,t)=0, n 2 = 1 (159) 
which is periodic in a\, h{a\ + 2n,t) = n(ai), i.e. 

SL(C,h) = lim L(C,C") . (160) 

In this last equation, 

is the familiar GauB linking number. 

The vector n(a\, t) defines a framing of the loop C(t), which implies a thick- 
ening of the loop into a ribbon bounded by C(t) and C'{t). Let C denote the 
loop along the middle of this ribbon which is generated by x(ai,t) + §n(ci, t) 
with t fixed. Taking the limit e —>■ 0, one finds for the self-linking number 



SL(C, h) = W r (C) + T W (C, h) (162) 

Here 

W r {C) = L{C,C) (163) 
is the writhing number and 

T w (C,n) = ida x [e{(J U t) x n{a x , t)\ ■ h\a u t) (164) 
Ztt J 
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denotes the twist, where 



e{a 1 ,t) = — r — (165) 

\x {<ri,t)\ 

is the unit vector tangential to the loop C(t). The twist T w {C,h) is not 
a topological invariant but depends on the choice of the framing n{a\,t). 
A global change of n(ai,t) changes T w (C,h) by some integer. The twist 
T w (C,h) expresses the torsion of the ribbon bounded by C(t) and C'(t) and 
represents the Polyakov spin factor in D = 3. It should be emphasized 
that, while the self-linking number and the twist depend on the choice of the 
framing 77.(01, t), their difference, i.e. the writhing number, is independent 
of this arbitrariness. Furthermore, it can be shown that the self-intersection 
number I2{S, S) is related to the difference in the self-linking number between 
the initial and final time t\ and ti of the four-dimensional space-time manifold 



h{S,S)= -SL{C{t),n{t))\^ ■ ( 166 ) 

In an analogous way, Polyakov's intersection number is related to the twist 
ofC(t), 

I 1 (S,S) = T w (C(t),n(t))\ t t= t 2 1 , (167) 
so that the total intersection number is given by the writhing number 

I(S,S)= -W r (C(t))\l=tl ■ (168) 

A transversal self-intersection point of the vortex sheet gives rise to a self- 
intersection point of the vortex string C(t) at some time t. If the loop C(t) 
self-intersects at some time t = i, the writhing number changes by 2 from 
the time before t until the time after the intersection. This contribution 
2 corresponds to the fact that each (transversal) self-intersection point of 
the vortex sheet contributes 2 units to the intersection number. Since for 
all closed vortex sheets, self-intersection points occur pairwise, the vortex 
string C(t) must self-intersect an even number of times, so that the SU(2) 
Pontryagin index 

1, 



u[A(S)\ = --W r (C(t)) 



t=t 2 



(169) 



t=ti 

becomes an integer, as observed already in section [| 
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C Topological charge of a thick SU(2) vortex 
intersection 

Consider two intersecting planar thick SU (2) vortices in four space-time di- 
mensions. By choice of coordinate system, the intersection region can be 
centered on the origin of space-time, and one of the vortices can be centered 
on the 1-2-plane, 

A« = T 3 fli f + f ] (0, 0, x 4 , -x 3 ) (170) 

cf. eq. ([H]), with an arbitrary profile function fi. On the other hand, let the 
other vortex be centered on the 3' — 4'-plane of a (primed) coordinate system 
related to the previous one by an arbitrary orthogonal transformation O, 

^' = T /^! + f\ x' 2 ,- X [A0) (171) 
x[ = OijXj (172) 

This also defines the transformation of the direction of the four-vector A^' 
when going to the unprimed coordinate system, = T A^'. In the 
unprimed coordinate system, therefore, the second thick vortex is given by 

,(2) _ rp f2[(OijXj) 2 + (02kX k ) 2 } (n n n n t \ ri7^ 

1 " 3 (0 llXl r + {0 2m x m y ^O^v-O^x,) (173) 

From this, one straightforwardly obtains the field strengths associated with 
the thick vortex fields (where the prime on the profile functions denotes 
differentiation w.r.t. the argument), 

= -F$ = d 3 A^-d,A^ = -2T 3 f[(x 2 3 + xl) (174) 

F$p = otherwise (175) 

and 

F^f = diAf-d^ (176) 

= 2T z {O l3 2l - 23 O ll )f 2 {{O lk x k ) 2 + {0 2l x l ) 2 ) ■ (177) 
The Pontryagin index is given by 

32^ / d4xe ^MFijF k i) (178) 
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v 



One can easily verify that the integrand vanishes for each thick vortex taken 
separately, 

e^MF^F^) = e^MHPH?) = (179) 

and therefore the only nonvanishing contribution to v results from the inter- 
section region of the two thick vortices, 

^- J *xti{F$F§>) (180) 



^ J d*x (0 12 21 - 22 O n ) f 2 ((O lk x k ) 2 + (0 2lXl ) 2 ) f[(x 



3 + X Aj 



2tt 2 

By the change of integration variables 

(x 1 ,x 2 ) — ► (x[, x' 2 ) = {0 lk x k , 2i xi) (181) 
this simplifies (up to the sign of the Jacobian) to 

H = — f dx[dx' 2 dx 3 dx 4 fl{x' 2 + x' 2 2 )f[{x 2 3 + x\) (182) 

27T J 

roo />oo 

= 2 dr dr'rr' ' f 2 (r' 2 ) f[(r 2 ) (183) 



o 



x 



= \f2(r'X- /i(O " (184) 

Due to the boundary conditions on profile functions (cf. the description after 
eq. (0)), one thus finally obtains 

M = \, (185) 

generalizing the result for thin 577(2) vortices derived in section |5|. The spe- 
cific form of the profile functions has completely canceled, and the present 
treatment thus in particular also includes the case of thin SU(2) vortices. 
Also all reference to the coordinate rotation O specifying the relative orien- 
tation of the vortex planes has vanished. 
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